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ABSTRACT 


The  research  program  during  the  contract  period 
consisted  of  the  analysis  of  the  CJbitron/FEL  amplifier  in 
three -dimensions .  The  principal  configuration  of  interest 
consisted  of  the  propagation  of  an  energetic  electron  beam 
through  a  loss  -free  rectangular  waveguide  in  the  presence  of 
a  linearly  polarized  wiggler  field  with  parabolically 
tapered  pole  pieces.  The  purpose  of  the  tapered  pole  faces 
is  to  provide  a  mechanism  for  focussing  the  electron  beam 
in  the  plane  of  the  bulk  wiggler  induced  oscillation.  A 
nonlinear  theory  and  simulation  code  has  been  developed  to 
study  this  configuration  which  can  treat  a  multiple  mode 
interaction,  harmonic  growth,  efficiency  enhancement  by 
means  of  a  tapered  wiggler,  the  effect  of  beam  thermal 
spread  on  the  interaction,  the  injection  of  the  beam  into 
the  wiggler,  and  detailed  facets  of  the  particle  dynamics 
such  as  Betatron  oscillations  and  velocity  shear. 
Comparisons  of  the  experiment  at  the  Lawrence  Livermore 
National  Laboratory  are  excellent.  Some  of  the  principal 
results  of  the  simulation  relate  to  the  harmonic  and  tapered 
wiggler  interactions.  In  the  case  of  the  harmonic 
interaction,  the  nonlinear  simulation  shows  both  high  gain 
and  good  efficiency  at  the  harmonics,  although  the 
interaction  is  increasingly  sensitive  to  the  effects  of 
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velocity  spread  at  the  higher  harmonics.  The  tapered 
wiggler  interaction  shows  high  efficiencies  are  possible 
and  that  one  effect  of  the  tapered  wiggler  is  to  reduce  the 
sensitivity  to  beam  velocity  spread.  An  additional  area  of 
study  is  the  enhancement  of  the  efficiency  of  a 
configuration  based  upon  an  electromagnetic  wave  wiggler  by 
means  of  the  tapering  of  an  external  axial  guide  field.  The 
purpose  of  the  electromagnetic  wave  wiggler  is  to  achieve  a 
higher  operating  frequency  for  a  given  beam  energy  (relative 
to  a  magnetostatic  wiggler) ,  however,  electromagnetic  wave 
wigglers  are  inherently  more  difficult  to  control  (i.e. 
taper)  for  efficiency  enhancement  purposes.  Thus,  the  use 
of  a  tapered  axial  guide  field  for  this  reason  may  have 
important  applications. 


1.  GENERAL  DISCUSSION 


The  fundamental  Interaction  mechanism  of  the  Ubitron  is 
the  same  as  that  of  the  Free-Electron  Laser  (FEL),  the 
distinction  between  the  two  is  that  the  term  Ubitron 
specifically  refers  to  a  microwave  tube  while  the  term  FEL 
is  often  used  to  describe  devices  operating  over  the  entire 
spectrum.  The  common  Ubitron/FEL  interaction  mechanism 
derives  from  the  propagation  of  an  energetic  (but  not 
necessary  relativistic)  electron  beam  through  a  rippled 
magnetic  field  (called  a  "wiggler"  or  "undulator" ) .  The 
effect  of  the  wiggler  field  is  to  induce  an  oscillation  in 
the  transverse  velocity  which  is  to  induce  an  oscillation  in 
the  transverse  velocity  which  gives  rise  to  spontaneous  or 
incoherent  radiation.  Coherent  radiation  arises  from 
electron  bunching  in  Uie  ponderomotive  potential  formed  by 
the  beating  of  the  radiation  and  wiggler  fields.  The 
wavelength  of  the  output  radiation  depends  on  both  the 
period  of  the  wiggler  field  and  the  energy  of  the  electron 
beam,  and  experiments  have  proven  the  mechanism  to  operate 
over  wavelengths  ranging  from  the  microwave  through  the 
visible  spectra.  Within  this  context,  therefore,  the 
Ubitron  represents  a  specific  application  of  the  concept  to 
the  microwave  regime  and  we  shall  use  the  term  to  describe 
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FEL's  which  employ  electron  beams  with  energies  less  than 
about  500  keV. 

The  primary  configuration  studied  during  the  contract 
period  consisted  of  the  propagation  of  an  electron  beam 
through  a  loss -free  rectangular  waveguide  in  the  presence  of 
a  linearly  polarized  wiggler  field.  The  principal 
motivation  for  this  work  derives  from  the  interest  in  Code 
6840  at  the  Naval  Research  Laboratory  in  LJbitrons  based  upon 
extremely  short  period  planar  wigglers  and/or  higher 
harmonic  Ubitrons .  The  advantage  to  be  derived  from  such 
configurations  is  the  possiblity  of  generating  high 
frequency  radiation  using  low  energy  electron  beams.  The 
detailed  wiggler  field  model  we  employ  describes  a  planar 
wiggler  with  parabol ically  tapered  pole  faces.  This  type  of 
wiggler  field  provides  for  enhanced  focussing  in  the 
direction  of  the  bulk  wiggler -induced  motion,  and  was  first 
employed  experimentally  by  Phillips1.  An  analytic  model  for 
this  type  of  field  has  been  given  by  Scharlemann2 ,  who  found 
that  the  field  can  be  represented  in  the  form 
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where  By,  denotes  the  amplitude  and  kw  (  =  2tt/\Vi,  where  Xw  is 
the  period)  represents  the  wavenumber. 

The  analytical  technique  employed  is  based  upon 
experience  gained  in  the  extensive  prior  treatment  of 
Ubitron/FEL  configurations  based  upon  helical  wiggler 
fields3-5.  Since  we  are  interested  in  an  amplifier  model 
only  single -frequency  wave  propagation  is  considered,  which 
permits  Maxwell's  equations  to  be  averaged  over  a  wave 
period.  As  a  result,  the  fast  time  scale  phenomena  are 
eliminated  from  the  problem.  In  addition,  only  a  "beamlet" 
of  electrons  which  enter  the  interaction  region  within  one 
wave  period  needs  to  be  included  in  the  analysis.  These  two 
benefits  result  in  major  improvements  in  computational 
efficiency  in  comparison  with  full-scale  particle -in -cell 
simulation  codes.  The  boundary  conditions  imposed  by  the 
waveguide  wall  may  be  treated  by  the  expansion  of  the 
radiation  field  in  terms  of  the  vacuum  waveguide  modes.  In 
previous  work  with  helical  wiggler  configurations  a 
cylindrical  waveguide  was  assumed;  however,  in  treating  the 
planar  wiggler,  we  have  assumed  the  boundary  conditions  to 
be  determined  by  a  rectangular  waveguide. 
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While  our  primary  interest  is  in  the  use  of  short  period 
wigglers  and  higher  harmonic  interactions  which  permit  high 
frequency  operation  with  moderate  energy  electron  beams,  the 
only  operational  experiment  against  which  our  nonlinear 
theory  can  be  compared  is  the  ELF  experiment  at  Lawrence 
Livermore  National  Laboratory6'7.  The  ELF  experiment  makes 
use  of  a  3.5  MeV  electron  beam  and  a  planar  wiggler  field 
with  a  9.8  cm  period  and  an  amplitude  of  as  much  as  4  kG  to 
obtain  operation  at  a  frequency  of  35  GHz.  As  a 
consequence,  the  initial  development  and  testing  of  the 
simulation  code  were  performed  for  this  experiment. 

The  initial  formulation  described  a  single-mode 
analysis,  and  gave  results  in  reasonable  agreement  with  the 
experiment.  The  analytical  treatment  and  numerical  results 
have  been  published  in  Physical  Review  A  (see  Appendix  I). 
Resonant  interaction  was  found  with  the  TE0i,  TE21>  TM2 1 
modes.  The  saturated  powers  found  for  the  TE  modes  were  of 
the  order  of  200  MW  which  is  comparable  to  that  observed 
experimentally.  The  coupling  to  the  TM  mode  is  weaker  than 
for  either  of  the  TE  modes,  and  both  the  gain  and  saturated 
power  are  substantially  lower  than  for  the  TE01  and  TE2i, 
modes . 
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The  principal  difference  between  the  nonlinear  analysis 
developed  under  the  present  contract  and  other  treatments, 
such  as  described  in  the  FRED  code  in  use  at  the  Lawrence 
Livermore  National  Laboratory7,  is  in  the  treatment  of  the 
particle  dynamics.  The  particle  orbits  are  treated  by  means 
of  a  reduced  set  of  orbit  equations  obtained  by  an  average 
of  the  Lorentz  force  equations  over  a  wiggler  period  in  FRED 
and  virtually  every  other  nonlinear  treatment  of  the 
Free-Electron  Laser.  The  average  of  the  orbit  equations  is 
not  performed  in  the  present  analysis,  however,  and  the  full 
set  of  Lorentz  force  equations  (i.e.,  in  three  dimensions) 
is  integrated  for  each  particle  in  the  simulation.  As  a 
result,  the  particle  dynamics  are  treated  more  accurately  in 
the  present  analysis  than  in  simulation  codes  such  as  FRED. 
This  is  important  in  the  description  of  the  planar  wiggler 
interaction  because  the  magnitude  of  the  bulk  wiggler 
velocity  is  not  constant  for  this  configuration,  and  an 
oscillation  occurs  in  the  axial  velocity  with  a  period  of 
half  the  wiggler  wavelength.  This  rapid  oscillation  is  lost 
when  the  orbit  equations  are  averaged  over  a  wiggler  period, 
but  manifests  itself  in  an  oscillation  in  both  the  wave 
power  and  phase  at  a  wavelength  of  Xw/2 .  While  this  does 
not  affect  the  bulk  growth  of  the  wave,  it  does  introduce  an 
uncertainty  in  the  measurement  of  the  axial  evolution  of  the 
wave  power  and  phase  which  may  be  as  high  as  10  -  2056. 
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One  advantage  of  the  more  accurate  treatment  of  the 
particle  dynamics  employed  in  the  present  analysis  is  that 
the  injection  of  the  beam  into  the  wiggler  may  be  modelled. 
The  injection  process  is  crucial  to  the  efficient  operation 
of  the  Free-Electron  Laser  because  poor  matching  of  the  beam 
into  the  wiggler  can  result  in  a  large  axial  energy  spread 
which  degrades  the  interaction  efficiency.  In  order  to 
minimize  this  problem  most  experimental  configurations 
employ  an  adiabatic  entry  taper  on  the  wiggler  amplitude 
which,  in  most  cases,  exceeds  four  wiggler  periods  in 
length.  Simulation  of  the  injection  process  shows  that  the 
adiabatic  injection  process  is  effective  for  entry  taper 
regions  longer  than  approximately  4-5  wiggler  periods. 
Entry  taper  regions  shorter  than  this  are  found  to  result  in 
substantial  declines  in  the  operating  efficiency. 

A  second  advantage  that  derives  from  the  improved 
treatment  of  the  particle  dynamics  is  that  the  analysis  is 
applicable  to  the  study  of  harmonic  radiation.  The  harmonic 
interaction  holds  promise  for  the  production  of  high 
frequency  radiation  using  electron  beams  of  moderate  energy. 
This  is  an  important  future  application  of  the  analysis. 
However,  in  experiments  designed  to  operate  in  the  visible 
spectrum,  the  harmonic  radiation  would  occur  in  the 
ultraviolet  which  can  result  in  fatal  damage  to  the  optical 
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system  As  a  result,  the  efficiencies  and  growth  rates  of 
the  harmonic  interaction  can  have  important  implications  for 
a  wide  variety  of  Free -Electron  Laser  applications.  The 
nonlinear  simulation  indicates  that  the  harmonic  interaction 
in  a  planar  wiggler  can  be  significant  with  growth  rates 
comparable  to  that  found  at  the  fundamental.  The  saturation 
efficiency  is  found  to  decrease  most  sharply  between  the 
fundamental  and  the  first  harmonic,  but  shows  a  weaK 
variation  with  harmonic  number  thereafter.  Thus,  the  total 
power  carried  by  the  harmonics  can  be  a  significant  fraction 
of  that  found  at  the  fundamental.  This  is  a  particularly 
encouraging  result  from  the  standpoint  of  the  design  of  a 
higher  harmonic  device.  However,  the  simulation  also  shows 
that  the  harmonics  are  substantially  more  sensitive  to  the 
effect  of  an  axial  energy  spread  than  the  fundamental,  and 
more  precise  control  of  the  electron  beam  will  be  required 
to  perform  a  harmonic  experiment.  These  results  have  been 
published  in  Physical  Review  A  (see  Appendix  II). 

Another  important  configuration  employed  in  the 
Free-Electron  Laser  interaction  mechanism  is  the  enhancement 
of  the  saturation  efficiency  by  means  of  a  tapered  wiggler. 
The  physical  basis  for  the  efficiency  enhancement  can  be 
explained  by  first  noting  that  wave -particle  resonance  in 
the  Free-Electron  laser  occurs  when  w  »  (k  +  kw)  vz,  where 
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s  vz  denotes  the  axial  electron  velocity,  kw  is  the  wiggler 

wavenumber,  and  (cj,  k)  denote  the  frequency  and  wavenumber 
of  the  radiation.  As  the  beam  gives  up  energy  to  the  waves, 
the  axial  electron  velocity  decreases  and  the  particles  drop 
out  of  resonance  with  the  wave.  The  purpose  of  the  wiggler 
tapering  is  to  accelerate  the  beam  in  the  axial  directon  so 
as  to  maintain  the  resonance  over  an  extended  interaction 
length.  This  can  be  accomplished  by  means  of  a  decrease  in 
the  wiggler  amplitude  which  affects  the  partition  of  energy 
between  motion  in  the  perpendicular  and  parallel  directions 
by  decreasing  the  perpendicular  velocity  of  the  bulk  wiggler 
motion  while  increasing  the  axial  velocity.  The  tapered 
wiggler  interaction  mechanism  has  been  addressed  in  the 
nonlinear  analysis  (see  Appendix  I),  with  particular 
emphasis  placed  on  the  effects  of  an  axial  energy  spread  on 
the  maximum  achievable  saturation  efficiency.  The  principal 
result  is  that  the  tapered  wiggler  interaction  is 
significantly  less  sensitive  to  the  effect  of  the  energy 
spread  than  is  the  case  lor  a  uniform  wiggler.  This  is 
encouraging  from  the  standpoint  of  the  design  of  a  tapered 
wiggler  harmonic  experiment,  since  the  higher  sensitivity  to 
beam  energy  spread  encountered  for  the  harmonic  interaction 
may  be  compensated  for  by  the  tapered  wiggler.  This  will  be 
an  important  area  of  future  study. 
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Recognizing  that  a  single -mode  representation  is  not 
appropriate  for  an  overmoded  configuration  in  which  the  beam 
is  resonant  with  more  than  one  mode,  we  undertook  to  include 
multiple  modes  in  the  simulation.  The  procedure  for  the 
inclusion  of  multiple  modes  within  the  present  formulation 
involves  (1)  the  integration  of  particle  trajectories  in  the 
presence  of  the  aggregate  field,  and  (2)  the  computation  of 
the  coupling  of  each  mode  to  the  total  current  to  obtain  the 
growth  of  the  wave  mode.  It  should  be  remarked  that 
although  multiple  modes  can  be  treated  within  the 
framework  of  the  analysis,  all  the  modes  are  at  the  same 
frequency  in  order  to  perform  the  average  of  the  field 
equations  over  a  wave  period.  The  results  of  the  analysis 
have  been  submitted  for  publication  in  the  Physical  Review  A 
(Appendix  III).  The  analysis  has  also  been  applied  to  the 
case  of  the  experiment  at  LLNL7  with  excellent  results.  The 
experiment  made  use  of  a  3.5  MeV/850  A  electron  beam  with 
an  initial  radius  of  1  cm  propagating  through  a  waveguide  of 
9.8  cm  x  2.9  cm  cross  section.  The  wiggler  field  was 
characterized  by  an  amplitude  of  3.72  kG  and  a  period  of  9.8 
cm.  The  input  radiation  source  was  a  magnetron  producing  50 
kW  at  a  frequency  of  34.6  GHz  predominantly  in  the  TE01  mode 
of  the  waveguide.  The  best  estimates  available8  indicate 
that,  perhaps,  1%  of  the  total  power  is  carried  by  the  TE2i 
and  TM2l  modes.  The  observed  output  power  for  a  uniform 
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wiggler  (i.e.  constant  amplitude)  was  180  MW,  and  rose  to 
1GW  (for  an  efficiency  of  34%)  when  the  amplitude  was 
decreased  by  55%  over  a  distance  of  1.1  m.  The  simulation 
code  for  these  parameters  yield  an  output  power  of  185  MW 
for  the  uniform  wiggler  case,  and  1  GW  for  the  tapered 
wiggler  interaction.  It  is  important  to  remark.,  however, 
that  the  theory  is  not  in  complete  correspondence  with  the 
experiment.  The  most  important  difference  lies  in  the  use 
of  tapered  pole  pieces  to  provide  for  beam  focussing  in  the 
plane  of  the  bulk  wiggler  motion  in  the  simulation,  while 
the  experiment  employed  a  quadrupole  magnetic  field  for  this 
purpose.  For  this  reason,  the  excellent  agreement  between 
the  theory  and  the  experiment  should  not  be  overstated. 
However,  it  does  indicate  that  the  essential  physics  of  the 
interaction  is  included  in  the  theory,  and  the  simulation 
will  now  be  turned  to  the  study  of  harmonic  interactions  in 
earnest . 

Contributed  papers  on  this  work  have  been  presented  at 
the  9th  Free  Electron  Laser  Conference  in  Williamsburg, 
Virginia  (9-13  September,  1987),  and  at  the  American 
Physical  Society/Division  of  Plasma  Physics  meeting  in  San 
Diego,  CA  (2-6  November,  1987).  In  addition,  recognition  of 
the  importance  of  these  contributions  to  the  understanding 
of  the  physics  of  the  Ubitron/Free  Electron  Laser  has 
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resulted  in  an  invited  paper  presented  at  the  9th 
International  Free -Electron  Laser  Conference. 

An  additional  area  of  study  is  the  configuration  based 
upon  electromagnetic  wave  wigglers.  The  advantage  of  this 
type  of  configuration  is  that  a  higher  frequency  interaction 
may  be  obtained  for  a  given  electron  beam  energy  relative  to 
a  magnetostatic  wiggler.  However,  one  disadvantage  of  the 
electromagnetic  wave  wiggler  is  that  it  is  more  difficult  to 
control  (than  magnetostatic  wigglers)  for  the  purpose  of 
efficiency  enhancement  through  a  tapered  amplitude  or 
period.  In  order  to  circumvent  this  problem,  we  have 
proposed  the  use  of  an  auxiliary  axial  guide  field.  The 
axial  guide  field  is  useful  from  the  standpoint  of  confining 
(or  focussing)  the  beam.  However,  the  nonlinear  interaction 
efficiency  may  be  enhanced  by  tapering  the  guide  field  as 
easily  as  tapering  the  wiggler  field4.  As  a  consequence,  we 
have  analyzed  the  effect  of  a  tapered  axial  guide  field  on 
the  electromagnetic  wave  wiggler  configuration,  and  the 
results  have  been  published  in  the  IEEE  Journal  of  Quantum 
Electronics  (see  Appendix  IV).  The  results  indicate  that 
the  tapered  axial  guide  field  may  have  important 
applications  to  the  enhancement  of  the  efficiency  in 
electromagnetic  wave  wiggler  configurations. 
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The  nonlinear  evolution  of  the  free-electron-laser  (FEL)  amplifier  is  investigated  numerically  for  a 
configuration  consisting  of  a  planar  wiggler  with  parabolically  tapered  pole  pieces.  A  set  of  coupled 
nonlinear  differential  equations  is  derived  in  three  dimensions  which  governs  the  self-consistent  evolu¬ 
tion  of  the  TE  and  TM  modes  in  a  loss-free  rectangular  waveguide  as  well  as  the  trajectories  of  an 
ensemble  of  electrons.  The  initial  conditions  are  chosen  to  model  the  injection  of  a  cylindrically  sym¬ 
metric  electron  beam  into  the  wiggler  by  means  of  a  region  with  an  adiabatically  tapered  wiggler  am¬ 
plitude,  and  the  effect  of  an  initial  beam  momentum  spread  is  included  in  the  formulation.  Both 
self-field  and  space-charge  effects  have  been  neglected,  and  the  analysis  is  valid  for  the  high-gain 
Compton  regime.  In  addition,  the  phase  stability  of  the  FEL  amplifier  against  fluctuations  in  the 
beam  voltage,  the  enhancement  of  the  efficiency  by  means  of  a  tapered  wiggler  amplitude,  and  har¬ 
monic  generation  are  also  studied.  Numerical  simulations  are  conducted  to  model  a  35-GHz 
amplifier  with  an  electron  beam  energy  of  3.5  MeV,  and  good  agreement  is  found  between  the  simu¬ 
lation  and  an  experiment  conducted  by  Orzechowski  and  co-workers.  Significantly,  the  results  indi¬ 
cate  that  a  tapered  wiggler  configuration  is  somewhat  less  sensitive  to  the  beam  thermal  spread  than  a 
uniform  wiggler  system. 


I.  INTRODUCTION 

The  free-electron  laser  (FEL)  and  the  ubitron1  have 
been  successfully  demonstrated  as  radiation  sources  over  a 
broad  frequency  range  from  the  microwave1  13  through 
the  optical14- 21  spectra.  The  distinction  between  the  ubi¬ 
tron  and  the  FEL  is  not  well  defined  in  the  literature; 
however,  we  find  it  convenient  to  differentiate  between  the 
FEL  and  the  ubitron  primarily  on  the  basis  of  the 
electron-beam  energy.  As  such,  we  refer  to  those  devices 
as  ubitrons  when  the  beam  energy  is  below  500  keV.  Al¬ 
though  this  definition  is  somewhat  arbitrary,  operation  at 
these  energies  generally  involves  frequencies  close  to  the 
waveguide  cutoff  of  the  device,  and  the  ubitron  may  be 
thought  of  as  a  weakly  relativistic  FEL  operated  as  a  mi¬ 
crowave  tube.  In  either  case,  however,  the  physical- 
interaction  mechanism  is  the  same,  and  relies  on  a  period¬ 
ically  rippled  magnetic  field  (referred  to  as  the  wiggler 
field)  to  induce  an  oscillatory  motion  in  the  electron 
beam.  The  interaction  between  the  transverse  component 
of  the  oscillatory  motion  and  the  radiation  field  results  in 
an  axial  bunching  of  the  electron  beam  which  is  the 
source  of  the  instability.  This  axial-bunching  mechanism 
can  be  thought  of  as  the  result  of  the  ponderomotive  po¬ 
tential  formed  by  the  beating  of  the  wiggler  and  radiation 
fields.  The  precise  form  for  the  wiggler  field  can  take  a 
variety  of  configurations,  and  ubitrons  and  FEL's  have 


been  constructed  using  both  helically2-810- 121415  and 
linearly'  9 13,16-21  polarized  wiggler  fields.  In  addition,  a 
wiggler  configuration  has  been  proposed  which  makes  use 
of  a  rotating  quadrupole  design.22 

The  motivation  for  the  present  work  is  to  develop  a 
nonlinear  theory  and  simulation  code  for  a  ubitron  or 
FEL  amplifier  based  upon  a  linearly  polarized  wiggler 
field.  The  advantage  of  a  planar  wiggler  design  over  that 
of  a  helical  configuration  is  ease  of  construction  and 
modification.  Linearly  polarized  wigglers  are  readily  con¬ 
structed  from  permanent  magnet  arrays  which  can  be 
easily  adjusted  to  a  tapered  wiggler  configuration.  In  this 
paper  we  derive  a  fully  three-dimensional  nonlinear 
analysis  of  the  FEL  and  the  ubitron  for  the  planar  wiggler 
configuration.  The  analysis  follows  that  described  previ¬ 
ously  for  a  helical  wiggler  configuration,23-27  and  involves 
the  derivation  of  a  set  of  coupled  nonlinear  differential 
equations  which  self-consistently  describe  the  evolution  of 
both  an  ensemble  of  electrons  and  the  electromagnetic 
fields  in  a  rectangular  waveguide.  Space-charge  fields  are 
neglected  in  the  analysis;  therefore,  the  treatment  is  ap¬ 
plicable  to  the  high-gain  Compton  (or  strong-pump)  re¬ 
gime.  The  nonlinear  current  which  mediates  the  interac¬ 
tion  is  computed  from  the  microscopic  behavior  of  the 
electrons  by  means  of  an  average  of  the  electron  phases 
relative  to  the  ponderomotive  wave  formed  by  the  beating 
of  the  wiggler  and  radiation  fields.  The  detailed  wiggler 
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model  we  employ  includes  the  effect  of  parabolically  ta¬ 
pered  pole  pieces  in  order  to  provide  for  electron  focusing 
in  the  plane  of  the  bulk  wiggler  motion,  Further,  the  in¬ 
jection  of  the  electron  beam  into  the  wiggler  held  is 
modeled  by  allowing  the  wiggler  amplitude  to  increase 
adiabatically  from  zero  to  a  constant  level.  The  procedure 
used  also  permits  the  inclusion  of  an  arbitrary  taper  of  the 
wiggler  amplitude  for  the  purpose  of  the  enhancement  of 
the  interaction  efficiency.  The  overlap  between  the  elec¬ 
tron  beam  and  the  transverse  mode  structure  of  the  TE 
and  TM  modes  is  included  in  a  seif-consistent  manner,  so 
that  no  arbitrary  ''filling  factor"  need  be  included  in  the 
analysis.  Since  the  problem  of  interest  is  that  of  a  FEL  or 
ubitron  amplifier,  only  single-frequency  propagation  is 
considered.  This  permits  an  average  over  a  wave  period  to 
be  performed  which  eliminates  the  fast-time-scale  phe¬ 
nomena  from  the  formulation,  and  results  in  a  great  in¬ 
crease  in  computational  efficiency  over  a  fuli-scai* 
particle-in-cell  simulation, 

In  organization  of  this  paper  is  as  follows.  The  general 
equations  are  derived  in  Sec.  II.  The  numerical  solution 
to  the  dynamical  equations  is  given  in  Sec.  IV,  in  which  a 
particular  example  is  treated  in  depth  which  corresponds 
to  a  recent  expenmen:  by  Orzechowski  and  co-workers. 1 3 
Three  distinct  waveguide  modes  are  found  to  grow  in 
simulation,  and  we  consider  each  of  these  modes  in  detail 
including  the  bandwidths  and  reiative  growth  rates  and 
saturation  efficiencies.  Also  considered  are  the  effects  of 


'll  venations  in  the  length  of  the  entry  region  on  the  satu¬ 
ration  efficiency,  (2;  an  initial  momentum  spread  in  the 
electron  beam  U.e.,  prior  to  the  injection  into  the  wigglen, 
(3)  the  scaling  of  the  efficiency  with  beam  current,  and  i4) 
the  enhancement  of  the  interaction  efficiency  through  a 
linearly  tapered  wiggler  amplitude,  These  issues  are  treat¬ 
ed  in  an  abstract  manner  in  Sec.  Ill  because  the  analysis 
does  not  correspond  to  all  aspects  of  the  experimental 
configuration.  However,  a  comparison  between  the  simu¬ 
lation  and  the  experiment  is  made  in  Sec.  IV  in  which  we 
give  a  summary  and  discussion. 


II.  GENERAL  FORMULATION 

The  configuration  we  consider  is  that  of  a  relativistic 
electron  beam  propagating  through  a  loss-free  rectangu¬ 
lar  waveguide  in  the  presence  of  a  linearly  polarized 
wiggler  magnetic  field.  The  wiggler-field  model  that  we 
employ  is  that  encountered  when  the  individual  magnets 
in  the  wiggler  have  parabolicallv  tapered  pole  faces, 
which  provides  for  electron-beam  focusing  in  the  plane 
of  the  principal  wiggler  motion.  This  technique  was  first 
employed  experimentally  by  Phillips.1  A  detai'ed 
analysis  of  the  magnetic  field  produced  by  a  wiggler  with 
parabolicaHy  tapered  pole  pieces  was  undertaken  by 
Scharlemann,;s  who  showed  that  the  wiggler  field  is  of 
the  form 


f  ,  f  .  f  kwx  i  ,  Ky  L 

B Jxi-IF.  i  cos (k.z)  ,  sinh  — ,  smh  —j=-  je. 


cosh  I  -7=-  cosh  '  — ■»-  ;  t, 

•  vi  i  v2  ; 


-t/2cosh  |~7T  i  s‘nh  !"7f  !  sinlA^zle,  f  . 

1  V1  J  l V2  J  i 


where  Bw  denotes  the  wiggler  amplitude,  and 
k.v.i  ~2tt/aw  )  is  the  wiggjer  wave  number.  We  model  the 
injection  of  the  electron  beam  into  the  wiggler  by  allowing 
the  wiggler  amplitude  to  increase  adiabatically  from  zero 
to  a  constant  level  over  -Vw  wiggler  periods  In  addition, 
since  we  intend  to  study  efficiency  enhancement  by  means 
of  a  tapered  wiggler,  the  wiggler  amplitude  will  be  tapered 
starting  at  some  point  z0  downstream  from  the  entry  re¬ 
gion  in  a  linear  fashion.  To  this  end,  we  choose 

|  Bw  sin^i  kw  z  /ASW  i,  0  <  z  <  Sm  t.w 
Bjz'.si  <£rk,  ywku.<2<,2i  d) 

\Bw[l +ewkw''Z -Zq>].  2>*o 


describes  the  slope  of  the  taper.  Since  the  fringing  fields 
associated  with  the  tapered  wiggler  amplitude  are  neglect¬ 
ed,  this  representation  is  strictly  valid  only  if  the  slope  of 
the  taper  is  small.  Within  the  entry  region  this  implies 
that  .V,,  must  be  large,  while  for  z  >  z0  w«  mu*t  require 


that  i  ew  ;  «  1 . 

The  boundary  conditions  at  the  waveguide  wall  may  be 
satisfied  by  expanding  the  vector  potential  in  terms  of  the 
orthogonal  basis  functions  of  the  vacuum  waveguide. 
Thus,  we  write  the  vector  potential  of  the  radiation  in  the 
form 

6 Aix,t}=  J'  5-4!Ku)ei,rU,.v)cosa  i4' 

n  =  0 

for  the  TE  modes,  and 


for  the  TM  modes,  where  for  frequency  m  and  wave  num¬ 
ber  k  (z), 


as  [‘ dz’klz')—a)t  . 

■'o 
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In  addition,  2’  indicates  that  /  and  n  are  not  both  zero, 
and 
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are  the  polarization  vectors.  In  this  representation  the 
waveguide  is  assumed  to  be  centered  at  the  origin  and 
bounded  by  —  a  /2  <  x  <a/2  and  —  b/2<y<b/2.  Asa 
consequence,  X  =x  +  a  /2,  Y=y+b/ 2,  and 

(,2  „2  1/2 

K  =  rr  -4  +  77  (9) 

a  b 

denotes  the  cutoff  wave  vector.  It  is  implicitly  assumed 


that  both  S-4|n(z)  and  k(z)  vary  slowly  over  a  wave 
period. 

The  microscopic  source  current  can  be  written  as  the 
following  sum  over  individual  particle  trajectories: 

,vr 

8J(x,t)=-e  2  v,(z;p,o,x,o,y,o,/,0) 

1  =  1 

X  5[  xt -  xu  (z ;  p,o,2C,0,y,o,  t,0 )  ] 

8[t  r,(z  ,Pio,x,o,yio,t,o)] 

X— — ; - —  ,  (10) 

where  L  is  the  length  of  the  interaction  region,  .Vr  is  the 
total  number  of  electrons,  nb  is  the  average  electron  densi¬ 
ty,  v,(z;p,0,jc,o,y,o,r,o)  is  the  velocity  of  the  ith  electron  at 
position  z  which  entered  the  interaction  region  (i.e., 
crossed  the  z  =0  plane)  at  time  tl0  and  transverse  position 
U,o,>-,o)  with  momentum  pl0,  and 

r, /;  ■  "" 

This  discrete  sum  over  particles  can  be  replaced  by  an  in¬ 
tegration  over  the  initial  conditions,  and  we  write 


8J (x,t)=—enb  f  dp0v20F0lp)  f  f  dx0dy0a  L(x0,y0)  f  ^  dt0a.l(t0)v(z;p0,x0,y0,t0) 

,  6[t-r(.r;p0,2c0.yo.'o)] 


X8[x1-x1(z;p0,jc0lj'0,to)]- 


where  uz0  is  the  initial  axial  velocity,  Ag  =  ab  is  the  area  of 
the  waveguide,  T  =  L/(v. 0),  and  a  L(x0,y0),  a  (fol>  ar>d 
Fn  (po>  describe  the  distributions  of  the  initial  conditions 
subject  to  the  normalizations  ( Ab  is  the  cross  sectional 
area  of  the  beam) 

f  f  dx0dynajx0,y0>=  Ab  ,  (13) 

fT/2  dt0a  (f0)=r  ,  (14) 

J  -  772 

and 

f  dp0F0(p0l  =  1  .  (15) 

Substitution  of  the  microscopic  fields  and  the  source 
current  into  Maxwell’s  equations  yields  the  equations 
which  govern  the  evolution  of  the  radiation  amplitude  and 
wave  number.  The  procedure  is  formally  identical  to  that 
described  previously  for  the  helical  wiggler,  cylindrical 
waveguide  configuration,23  and  involves  a  modal  orthogo- 
nalization  in  the  transverse  coordinates.  In  addition,  a 
quasistatic  assumption  is  made  in  the  sense  that  particles 
which  enter  the  interaction  region  at  times  t0  separated  by 
integral  multiples  of  a  wave  period  are  assumed  t.>  execute 
identical  trajectories.  As  a  result,  viz  ;po,x0.>'o,to 
+  2trAf/<u)  =  v(z;p0,jc0,>'o,fo)  for  ‘he  integer  .V,  and  a 
time  average  over  a  wave  period  can  be  performed  which 


permits  consideration  of  a  single  “beamlet”  that  includes 
electrons  which  enter  the  interaction  region  within  one 
ponderomotive  (or  wave)  period.29  This  greatly  improves 
the  computational  efficiency.  It  should  be  remarked  that 
the  ponderomotive  phase  of  each  particle  is  followed  self- 
consistently  in  the  analysis,  so  that  while  the  beamlet  is 
initially  distributed  over  one  ponderomotive  wave,  the 
final  state  may  describe  an  electron  distribution  which  has 
become  trapped  within  several  ponderomotive  buckets.  In 
view  of  this,  the  equations  which  govern  the  evolution  of 
the  TE|n  mode  are 

yT8crln+  ~-k2-kfn  5 o,n 
dz  c 


2k,/2-y-(k  l/26aln)  = 
dz 


_  /  COSO L  m  \ 

~ 8  c 2  ^ln\  1  1  7’ 

/  sing  ,„  \ 

~8c2  ln\  Me,n  7’ 


where  ba\n=ebA\„/mc2,  a)l=4ire2nb/m,  v  is  the  instan¬ 
taneous  velocity,  and  F|„  = }  when  either  1=0  or  n  =0 
and  unity  otherwise.  For  the  TM)n  mode  we  obtain  a 
similar  result, 
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Observe  that  theie  is  no  nontrivial  TM  mode  solution  when  either  /~ 0  or  n  ^ 0 . 

The  averaging  operator  <(  •  ■;>  is  defined  over  the  initial  conditions  of ‘he  beam  -note  that  the  instantaneous  posi¬ 
tions  and  momenta  of  the  electrons  are  implicit  functions  of  the  initial  conditions'1  and  includes  toe  effect  of  an  initial 
momentum  spread  by  means  of  the  distribution  function 


F-JpCi^A  exp[  - lp,c-po  'VA/>;]8!/>5 -p f0  -p;0 Wp* ■  , 


<20i 


where  p0  and  \p.  descriae  the  itutiai  bulk  momentum  and  momentum  spread  of  the  ream,  fj'x'  is  the  Heaviside  fun: 
non  and  the  normalization  constant  is 


A~  ~  f*°  dp,0e*p[-ipa-Pz)l/Zp;]  ' 


2 1 : 


Observe  that  the  distribution  is  monoenergetie,  but  contains  a  pitch  angle  spread  which  describes  an  axial  energy  spread 
■as  well  as  a  transverse  energy  spread!  given  approximately  by 


Am 
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1  -t-  21  —  1 


1 22/ 


P  o 


where  yc~‘  1  — p;  /m 2c‘:;  There  is  no  fundamental  difficulty  in  the  inclusion  of  an  overall  energy  spread  in  the 
analysis;  however,  the  additional  degree  of  freedom  requires  an  increased  number  of  panicles  in  the  Simulation.  The  6 
function  allows  us  to  perform  one  of  the  momentum-space  integrals  analytically,  and  to  write  the  averaging  operator  m 
the  form 


i>5-~£  /;’<Wo  /*5 4Psofto«*Pl  -'>j0  -po f  -Vo I  /:;>x0  f*'^dy0ojxCly 


iff  ■' 


(23) 


where  d'o !s— w»v  is  the  initial  ponderomotive  phase, 
60S  tM"::pyo/Pxih  and  0jO-=isO/c.  It  is  important  to 
recognize  that  this  average  includes  the  effect  of  the  over¬ 
lap  between  the  electron  beam  and  the  transverse  mode 
structure  of  the  radiation  field  in  a  seif-consistem  way. 

The  phase  variation  of  each  mode  can  be  analyzed  oy 
the  addition  of  an  equation  to  integrate  the  relative  phase 

<t>(*j“  f* dz'{k(z')-kQ]  ,  124) 

where  A,0S!W:/£’2-/c!Jn;1'J  is  the  wave  number  of  the 
vacuum  guide  Since  the  departure  of  k  (z)  from  the  vacu¬ 
um  wave  number  describes  the  effect  of  the  wave-particle 
interaction.  represents  a  measure  of  the  dielectric 
effect  of  the  FEL  interaction.  Thus,  we  integrate  the  ad¬ 
ditional  equation 

>-k-*0,  i23) 

dz 

for  both  the  TE  and  TM  modes. 

In  order  to  complete  the  formulation,  the  electron-orbit 
equations  must  also  be  specified.  Since  we  describe  an 
amplifier  model,  we  choose  to  integrate  in  1  and  write  the 
Lorentz  force  equations  in  the  form 

u,  -—p^  —  e8Ein- “VX'BU ,  (26) 

dz  c 


where  IF  is  given  bv  Eqs.  il)  and  -Z)  and  the  radiation 
fields  are  given  by  the  vector  potentials 

1  6  B 


Finally,  the  electron  coordinates  obey  the  equations 
d 

V‘TzX^' 


V‘dz>'~V-' 


and 

dii  1  w 

6  —  k  -i-  kw  — 

dz  Uj 

describes  the  evolution  of  ponderomotive  phase 


if ~ ibi] 4-  fldz' 


k  — 


0) 


.27; 

(23: 

29) 

(301 

(31) 


III.  NUMERICAL  ANALYSIS 

The  set  of  coupled  differential  equations  described  in 
Sec.  II  is  solved  for  an  amplifier  configuration  in  which  a 
single  mode  of  frequency  o  is  injected  into  the  system  at 
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z=0.  The  dynamical  equations  [Eqs.  ( 16)— ( 19)  and  (25)] 
for  the  fields  can  be  reduced  to  a  set  of  four  first-order 
ordinary  differential  equations  for  5<J|n, 
r^fs^r'ddnSaiJ/cfe],  k ,  and  the  relative  phase  4*. 
Hence,  the  numerical  resolution  of  the  problem  consists  in 
the  simultaneous  solution  of  6^7-  +  4  first-order  ordinary 
differential  equations,  where  .Vr  is  the  number  of  elec¬ 
trons.  The  algorithm  we  employ  is  the  fourth-order 
Runge-Kutta-Gill  technique.  While  this  technique  is 
somewhat  less  stable  than  the  fourth-order  Adams- 
Moulton  predictor-corrector  scheme,  it  has  the  advantage 
of  being  less  memory  intensive.  Indeed,  the  requirements 
placed  on  the  available  size  of  computer  memories 
represents  a  critical  practical  limitation  when  momentum 
spread  is  included.  The  averages  in  Eqs.  (16)-tl9)  are 
performed  by  means  of  an  jVth-order  Gaussian  quadrature 
technique  in  each  of  the  variables  (.'c0.>'n.iA>P--o-<*o)- 
The  initial  conditions  on  the  radiation  field  are  chosen 
such  that  r,n(2=0)=0,  k(z=0)  =  kn,  and  <t><2=0)  =  0 
for  an  arbitrary  initial  power  level.  Observe  that  the 
time-averaged  Poynting  flux  P[n  for  the  waveguide  modes 
is  related  to  the  field  amplitude  by  the  relation 


„  m'c  ab  ,  Q  2 
^0  =  7— r— r-wfcoa 
32e‘  ~r,n 

for  the  TE,n  mode,  and 
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for  the  TM|n  mode.  The  initial  state  of  the  electron  beam 
is  chosen  to  model  the  injection  of  a  continuous,  axisym- 
metric  electron  beam  with  a  uniform  cross  section;  hence, 
we  choose  <7  =  I  for  —ir<i'0<~  and  0  for  r0  <  Rh . 

The  particular  example  we  analyze  is  that  of  a  35-GHz 
amplifier  employing  an  electron  beam  with  an  energy  of 
3.5  MeV,  a  current  of  800  A.  and  an  initial  radius  of  1.0 
cm  which  propagates  through  a  waveguide  characterized 
by  a  =9.8  cm  and  b—  2.9  cm.  In  order  to  obtain  peak 
growth  rates  in  the  vicinity  of  35  GHz  we  choose  a 
wiggler  field  with  an  amplitude  of  3.72  kG  and  a  period  of 
9.8  cm.  with  an  entry  taper  of  .V,  =  5.  For  purposes  of  il¬ 
lustration,  the  first  case  we  consider  is  that  of  a  beam  with 
zero  momentum  spread  (A =0).  For  all  cases  discussed 
in  this  work  with  \y!  =0  the  choice  of  a  tenth-order 
Gaussian  algorithm  in  each  of  the  coordinates  (1 l’n,r0,d0) 
was  found  to  provide  an  accuracy  of  the  order  of  0.1%. 
The  initial  electron  distributions  in  the  axial  phase  space 
and  beam  cross  section  are  shown  in  Fig.  1 .  Each  dot  in 
the  illustration  of  the  axial  phase  space  describes  a  phase 
v heel  composed  of  100  electrons  distributed  throughout 
the  cross  section  of  the  electron  beam.  Each  phase  sheet. 
therefore,  represents  a  cross-sectional  slice  of  the  electron 
beam,  which  is  chosen  initially  as  shown  in  Fig.  Ilbt.  Al¬ 
though  each  phase  sheet  is  initially  chosen  to  be  identical, 
the  subsequent  evolution  of  the  electron  trajectories  in  the 
presence  of  the  radiation  field  is  followed  self-consistently. 
It  should  also  be  remarked  that  the  distribution  shown 
represents  a  uniform  electron  beam.  The  nonuniformity 
in  the  positions  of  the  electrons  chosen  by  means  of  the 


Axial  Phase  Space  (kwz  =  0) 

- - 1 - - — 1 - - 

(a) 


-rr  0  7 

w 


Beam  Cross  Section  (kwz  =  0) 

b  _  _  _ 

(b) 


y  0-- 


b  _  _ 

0  12a  0  0  12a 

X 

FIG.  1.  Initialization  of  1  at  the  axial  phase  space,  and  > bi  the 
beam  cross  section.  Each  point  in  the  axial  phase  space 
represents  the  superposition  of  100  particles  distributed 
throughout  the  cross  section  of  the  beam. 


Gaussian  algorithm  is  compensated  for  by  a  nonuniformi¬ 
ty  in  the  weighting  of  the  electrons. 

There  are  three  modes  which  are  resonant  in  th<*  vi¬ 
cinity  of  35  GHz;  specifically,  the  TE0|,  TE2,,  and  TM2, 
modes.  The  analysis  presented  in  this  work  deals  with 
wave-particle  interactions  with  single  modes,  and  we 
shall  deal  with  each  of  these  modes  individually.  The 
first  mode  we  treat  is  the  TE0|  mode  which  is  the 
lowest-order  resonant  mode  at  the  frequency  of  interest. 
The  detailed  evolution  of  the  wave  power  as  a  function 
of  axial  position  is  shown  in  Fig.  2  for  the  injection  of  a 
50-kW  signal  at  a  frequency  of  a i/cku.  =  11.3  (34.6  GHzi. 
As  shown  in  the  figure,  saturation  occurs  at  ku.z -  115 
1 1.79  m)  at  a  power  level  of  approximately  214  MW 
which  corresponds  to  an  efficiency  of  27  ^7. 75%.  Wave 
amplification  occurs  principally  within  the  uniform 
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wiggier  region  <z  >  :  which  is  1.30  m  for  this  case 

arid  yields  an  average  gain  of  about  2s  dB/m.  Although 
this  value  for  the  average  gain  is  lower  than  that  ob¬ 
served  in  some  other  experiments  (a  gain  of  approxi¬ 
mately  120  dB/m  has  been  observed  by  Gold  and  co¬ 
workers3!,  the  average  normaitzec  growth  rate 
IV*,,  ^0.05  is  quite  high  and  is  attributable  to  the  rel¬ 
atively  high  wiggier  amplitude  and  long  wiggier  period. 

One  feature  of  interest  shown  in  Fig.  2  which  merits 
some  discussion  ;$  the  oscillation  in  the  instantaneous 
power  which  occurs  with  a  period  of  This  is  nor 

found  for  the  case  of  helical  wiggier  configurations  for 
which  the  bulk  transverse  wiggier  motion  describes  a  helix 
with  a  transverse  velocity  of  relatively  constant  magni¬ 
tude.  In  contrast,  a  linearly  po.arized  wiggier  will  induce 
3  bulk  wiggier  motion  in  the  plane  normal  to  that  of  the 
wiggier  field  characterized  oy  an  oscillatory’  velocity.  In 
order  to  illustrate  this  qualitatively,  we  observe  that  for 
the  present  configuration  the  bulk  transverse  wiggier 
motion  is  aligned  along  the  a  axis  and  varies  approximate¬ 
ly  as 

n„. 

v  -  — f —  sintfc  z  ST  .  .34' 

*  Y)K. 

The  source  terms  contained  ir  the  dynamical  equations 
'16:-:  ;9;  are  derived  essentially  from  a  calculation  of 
v  J-5E 3 ;  hence,  the  principal  wave-particle  coupling  is 
with  the  x  component  of  the  radiation  field.  If  we  assume 
that  &EX  ~b£t  iini/a—un,  then  it  is  evident  that 
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The  interaction  occurs  when  the  ponderomotive  phase  is  a 
slowly  varying  function  of  axial  position  [i.e.. 


. a. -k  -He 


],  and  the  upper  beat  wave  which  varies  as 


cosi>  describes  wave  air.pliflcatton,  The  lower  beat  wave 
describes  the  oscillation  at  half  the  wiggier  period,  Al¬ 
though  the  spatial  average  of  the  contribution  of  the  lower 
beat  wave  vanishes  and  this  term  provides  no  contribution 
to  the  bulk  growth  of  the  wave,  the  instantaneous  values 


of  the  power  are  affected.  Indeed,  the  instantaneous  varia¬ 
tion  oi  the  relative  pnase  also  exhibits  an  oscillation  at 
half  the  wiggier  period. 

A  full  spectrum  of  the  TE,V  mode  is  shown  in  Fig.  i  in 
which  we  plot  the  saturation  efficiency  and  distance  to 
saturation  as  a  function  of  frequency  within  the  unstable 
band.  As  shown  in  the  figure,  wave  amplification  is  found 
for  frequencies  extending  from  ra  /cfc*  a  10  (30.6  GHz! 
through  ic/c.k.,  ~  14.2  (43.5  GHz-  with  a  peak  efficiency 
of  the  order  of  approximately  9.&%.  The  peak  growth 
rate  .as  measured  by  the  distance  to  saturation)  occurs  for 
a/eHr*  » !  2. 3  .37.7  GHz;,  which  is  somewhat  higher  than 
the  targeted  35  GHz.  However,  the  gain  bandwidth  is 
sufficiently  broad  that  the  growth  rate  Has  not  decreased 
sigrificaim  v  from  the  peak  value. 

The  variation  of  the  relative  phase  versus  axial  position 
is  li.ustrated  in  Fig.  4  for  oj/ck±  •=  10.4,  10.7,  11.0,  11.3, 
and  11.9.  As  is  evident  from  the  figure,  the  oscillation  at 
one  half  the  wiggier  period  is  also  manifested  in  the  rela¬ 
tive  phase.  The  bulk  variation  Ye.,  averaged  over  a 
wiggier  period-  shows  the  same  qualitative  behavior  as 
that  found  for  a  helical  wiggier. 2°  Specifically,  for  fre¬ 
quencies  at  the  low  end  of  the  gain  band  the  relative  phase 
decreases  up  ro  a  point  just  short  of  the  position  at  which 
the  power  saturates  (indicated  in  the  figure  by  an  arrow,, 
after  which  the  phase  remains  relatively  constant.  As  the 
frequency  increases,  the  variation  in  the  relative  phase  de¬ 
creases  until  t  critical  frequency  is  reached  iw/ckk  all 
for  the  particular  case  under  consideration)  for  which  the 
phase  is  found  to  vary  little  over  the  course  of  the  interac¬ 
tion.  This  unheal  frequency  is  typically  found  to  be  about 
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FIO.  2.  Plot  of  the  growth  of  the  TE:,  mode  with  ax-ai  pom 
lion. 
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FIG  2  Graph  showing  the  distance  to  saturation  and  the 
efficiency  of  the  TE_  mode  as  a  function  of  frequency. 
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TE0i  Mode  (a  =  9.8  cm;  b  =  2.9  cm;  Pin  =  50  kW) 


FIG.  4.  Plot  of  the  relative  phase  of  the  TE„  mode  vs  axial 
position  for  a> /ck,,  —  10.4,  10.7,  11.0,  11.3,  and  11.9. 


\0%  below  the  frequency  of  peak  growth.  For  frequen¬ 
cies  above  the  critical  frequency,  the  bulk  phase  increases 
monotonically.  It  should  be  remarked  that  the  gain  band 
under  discussion  corresponds  to  the  upper  trequency  in¬ 
tersection  between  the  beam  resonance  line  [<o  =  (k 
-i-ku.  H'.  ]  and  the  vacuum  waveguide  dispersion  curve. 
There  is  also  a  gain  band  associated  with  the  lower- 
frequency  intersection  which  shows  a  similar  variation  in 
the  relative  phase,  except,  that  there  is  an  inversion  in  the 
frequency  dependence.'6 

The  question  of  the  effect  of  the  mjection  proccss  of  the 
electron  beam  can  be  addressed  by  varying  the  length  of 
the  entry  taper  region.  The  results  of  this  analysis  are 
shown  in  Fig.  5  in  which  we  plot  the  saturation  efficiency 
and  the  distance  to  saturation  as  functions  of  the  length  of 
the  entry  taper  region  for  ,Va  >  3.  We  have  arbitrarily 
chosen  the  minimum  length  of  the  entry  taper  region  to 
be  three  wiggler  periods  since  the  fringing  fields  associated 
with  the  tapered  wiggler  field  have  been  neglected,  and  we 
feel  that  below  this  value  the  fringing  fields  will  be  impor¬ 
tant.  The  results  indicate  that  the  saturation  efficiency  in¬ 
creases  markedly  as  the  length  of  the  entry  region  in¬ 
creases  from  3Aa.  to  approximately  6 A.,,.,  after  which  the 
increase  in  the  efficiency  becomes  more  gradual.  We  attri¬ 
bute  this  increase  in  the  efficiency  to  a  decrease  in  the 
effective  momentum  spread  induced  by  the  injection  of 
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FIG.  5.  Plot  of  the  distance  to  saturation  and  the  efficiency  of 
the  TEn  mode  as  a  function  of  the  length  of  the  entry  taper  re¬ 
gion  for  ui/ckw  =  11.3. 


the  beam  into  the  increasing  wiggler  field.  It  should  also 
be  noted  that  the  increase  in  the  distance  to  saturation  is 
roughly  linear  for  \u.  >  5  and  corresponds  to  the  increase 
in  the  length  of  the  entry  region  (i.e.,  the  length  of  the 
uniform  wiggler  region  remains  relatively  constant). 

The  particle  dynamics  during  the  course  of  the  TE;i| 
mode  interaction  are  illustrated  in  Figs.  6  and  7.  In  the 
first  place,  we  remark  that  saturation  proceeds  by  means 
of  the  phase  trapping  of  the  beam  in  the  ponderomotne 
potential.  This  is  shown  in  Fig.  6  in  which  we  plot  the 
phase-space  distribution  of  the  beam  at  saturation.  The 
dashed  line  in  the  figure  represents  an  approximate 
separatrix  calculated  for  particles  at  the  beam  center; 
hence,  many  of  the  electrons  which  appear  outside  of  the 
separatrix  may  instead  be  on  trapped  particle  orbits  at  the 
edge  of  the  beam.  The  cross-sectional  evolution  of  the 
beam  is  shown  in  Fig.  7.  The  cross-sectional  projection  of 
the  beam  at  the  start  of  the  uniform  wiggler  region  (i.e., 
A-U  z  =  31)  is  shown  in  Fig.  7(a).  The  bulk  motion  of  the 
beam  exhibits  four  essential  features.  The  first  is  the  pri¬ 
mary  wiggler-induced  oscillatory  notion  which  shifts  the 
center  of  the  beam  off  axis  in  the  x  direction,  and  this 
shift  is  clearly  shown  in  Fig.  7(a).  The  second  feature  is 
that  the  transverse  wiggler  gradient  introduces  a  betatron 
oscillation  which  causes  a  macroscopic  scalloping  of  the 
beam  envelope.  In  addition,  on  a  microscopic  level  the 
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individual  electrons  come  into  a  focus  and  out  again  on 
*he  opposite  side  of  the  beam.  This  becomes  evident  ;n 
Fig.  ■  ai  by  noticing  that  the  two  closely  spaced  "spokes" 
are  oriemed  in  the  negative-*  direction  while  in  the  initial 
state  [rig.  Ob']  these  same  spokes  were  onemeb  in  the 
positive-*  direction.  The  third  feature  is  that  the  trans- 
ve;se  wiggle:  gradient  also  has  a  focusing  effect  on  the 
beam  which  results  in  a  reduction  in  the  maximum  beam 
radius  relative  to  the  initial  state.  The  betatron  Oscillation 
occurs  over  a  length  of  approximately  Azff-3  OfO.._  i.e  . 
k  19.  J:,  and  the  evolution  of  the  beam  cross  section 

over  this  distance  is  shown  in  Fig.  7  'from  /cc  =  3i-5;  \ 
Lastly,  the  geometry  of  the  wjggier  and  the  transverse 
gradients  tends  tc  distort  the  beam  into  an  elliptical  cross 
section  While  this  is  evident  in  Fig.  7  at  an  early  stage  of 
the  wiggle:,  it  is  shown  even  more  dramatically  in  Fig 
7!fi  which  shows  the  cross  section  at  saturation 
•  A-,r  =■  1 15 

We  now  consider  the  TE-,  mode,  and  plot  the  evolu¬ 
tion  of  the  wave  power  versus  axial  position  in  Fig.  8  for  a 
50-kW  input  signal  at  u/cku,  —  .1.3  ana  an  r'ectron  beam 
with  a  aero  initial  momentum  spread.  As  is  evident  from 
the  figure,  the  power  saturates  at  kxz  ^  104  at  about  194 
MW  for  an  efficiency  tji-,6.85%  In  comparison  with  the 
TEci  mode,  therefore,  we  conclude  that  the  average 
growth  rate  is  somewhat  higher  and  the  efficiency  lower 
for  the  TE;,  mode  at  this  frequency  A  complete  spec¬ 
trum  for  the  TEji  mode  is  shown  in  Fig.  9  in  which  we 
plot  the  distance  to  saturation  and  the  efficiency  versus 
frequency.  As  shown  in  the  figure,  gam  is  found  for  fre¬ 
quencies  ranging  from  v/ck^.  a-  8.9  through  a)/ckx  £  14 
with  a  peak  efficiency  of  approximately  12%.  As  a  result, 
both  the  bandwidth  and  peak  efficiency  are  higher  for  the 
TE;,  mode  than  for  the  TEn,  mode.  In  addition,  while 
this  frequency  falls  slightly  below  the  frequency  of  peak 


growth  rate  for  the  TEj,  mode,  it  is  close  to  the  frequency 
of  maximum  growth  rate  for  the  IEj:  mode. 

The  evolution  of  the  wave  powei  versus  tuuai  position 
is  shown  in  Fig.  10  foi  a  50-kW  input  signal  in  the  TM;, 
model  at  u,  ck.e  — 11.3  mid  an  electron  beam  W'it'n  zero 
initial  momentum  spread.  As  shown  in  the  figure,  the 
power  saturates  at  kxz~22~  with  approximately  68.5 
MW  for  an  efficiency  tj^2.45%.  This  is  a  much  lower 
growth  rate  and  efficiency  than  found  for  either  the  7E,v 
or  TF.2[  modes,  despite  the  fact  that  the  cutoff  frequency 
and  dispersion  curves  are  degenerate  for  the  TM;,  and 
TE;;  modes  The  difference  between  the  two  modes  lies 
in  the  transverse  mode  structure.  As  mentioned  previ¬ 
ously,  the  principal  component  of  the  wiggler-induced 
motion  is  aligned  with  the  *  axis;  hence,  the  wave- 
particle  interaction  is  governed  largely  by  the  x  com¬ 
ponent  of  the  electric  field.  Comparison  of  the  mode 
structures  for  the  TE;i  and  TM;!  modes  given  by  Eqs. 
i4!  and  (3;  shows  that  for  a  giver,  mode  amplitude 
16.4;,  i  the  ratio  of  the  x  component  of  the  electric  field 
of  the  TM;,  ro  that  of  the  TE?1  mode  is  approximately 
26/e~0. 59.  As  a  result,  the  wave-particle  coupling  is 
weaker  for  the  TM2,  mode.  Note  that  this  conclusion 
would  be  reversed  if  the  ccm* — ation  were  altered  such 
that  the  principal  component  of  '.he  wtggler  motion  were 
aligned  with  th-  •.  axis.  However,  one  effect  arising  from 
the  degert  Acy  of  the  dispersion  curves  is  that  the  fre- 
queo.,  of  interest  ‘U'/ckw~ !  1. 3 ;  ties  near  to  peak 
grow  h  for  bo»F  TE ; ;  and  TM2,  modes  This  is  shown 
clearly  in  Fig.  li  m  which  we  plot  the  distance  to  satu¬ 
ration  and  the  saturation  efficiency  as  a  function  of  fre¬ 
quency  foi  the  TM;,  mode.  As  might  be  expected,  the 
reduced  wave-particle  coupling  and  growth  rates  for  the 
TM;i  mode  result  m  a  narrower  bandwidth  for  instabili¬ 
ty,  and  we  obtain  wave  growth  for  frequencies  ranging 
from,  ai /ckw  ^  1 0. 2  - 1 2. 

The  question  of  the  variation  in  the  pha^e  of  FEL 
amplifiers  has  important  implications  for  many  of  the 
potential  applications  of  these  ceviees.  Since  high 
powers  and  efficiencies  have  been  demonstrated  experi- 
menially,15  one  such  application  may  be  as  a  high-power 
microwave  source  for  the  next  generation  of  radio¬ 
frequency  electron  accelerators.  However,  one  require¬ 
ment  for  such  an  application  is  good  phase  stability  of 
the  output  radiation  against  fluctuations  in  the  electrcn- 
beam  volrage.  Examination  of  Fig.  4  shows  that  the  rel¬ 
ative  phase  at  saturation  varies  rapidly  with  frequency. 
Since  there  is  a  correspondence  between  variations  in  the 
frequency  at  fixed  voltage  and  variations  in  the  voltage 
at  fixed  frequency,  we  might  expect  the  phase  at  satura 
tion  to  vary  rapidiy  with  beam  voltage.  This  is  indeed 
the  case  as  shown  in  Fig.  12.  in  which  we  plot  the  rela¬ 
tive  phase  at  a  fixed  axial  position  (chosen  to  correspond 
to  the  saturation  point  for  a  beam  voltage  of  3.5  MeV> 
versus  beam  voltage  for  the  TE?!,  TE21,  and  TM2; 
modes.  The  variation  in  the  relative  phase  is  approxi¬ 
mately  5 1 V 1  %  variation  in  the  beam  voltage  for  the 
TE(;t  mode,  43  V  l  %  variation  tn  the  beam  voltage  for 
the  TE7I  mode,  and  89  V 1  %  variation  in  the  beam  volt¬ 
age  for  the  TM;|  mode.  Such  rapid  variations  in  the 
phase  arc  consistent  with  results  obtained  previously  for 
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FIG-  3.  Plot  of  the  growth  of  the  TE?i  mode  with  axi&l  posi¬ 
tion 
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FIG  10.  Plot  of  the  growth  of  the  TMJ:  mode  with  axial  po¬ 
sition. 


a  helical  wiggler,  cylindrical  waveguide  configuration,^ 
and  contrasts  with  8V1%  variation  m  beam  voltage  for 
a  typical  Stanford  Linear  Accelerator  tSLAC!  klystron.30 
We  conclude,  therefore,  that  applications  which  require 
an  extremely  phase-stable  microwave  source  will  aiso  re¬ 
quire  an  electron-beam  source  with  a  very  low  level  of 
voltage  fluctuations. 

The  effect  of  an  initial  momentum  spread  is  shown  in 
Fig.  13  in  which  we  pic:  the  efficiency  versus  Ay,  /yc,  Lt 
the  TE0,,  TEi),  and  TM:1  modes.  From  the  beam- 


resonance  condition,  it  is  apparent  that  the  transition  to  a 
thermally  dominated  regime  occurs  when  At1,  /  r0  —  T ■„  / 
( k  -t ‘ku ).  Making  use  of  Eq.  (22)  we  And  that  this  transi¬ 
tion  occurs  at  about  Ay,  /)'a==  I89fc  for  the  TEg,  and  TEj< 
modes,  and  Ay,  /yQ^.S%  for  the  TMi:  mode.  Asa  conse¬ 
quence,  the  cases  shown  in  Fig.  19  are  well  short  of  the 
thermal-beam  regime.  The  efficiency  is  found  to  decrease 
in  an  approximately  linear  fashion  with  increasing  Ay,  for 
each  of  these  modes,  which  corresponds  with  results  ob- 
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tained  for  a  helical  wiggler  configuration.25  Note,  however 
that  the  present  results  were  obtained  with  the  same  mod¬ 
el  distribution  (20)  as  for  the  helical  wiggler,  and  the  de¬ 
tailed  variation  in  the  efficiency  with  momentum  spread 
can  be  expected  to  vary  with  the  choice  of  distribution.  Be 
that  as  it  may,  we  find  that  the  efficiency  drops  from 
r?~7.75%  to  17 ~ 4. 98%  for  the  TEqi  mode  as  Ay2  /y0  in¬ 
creases  to  about  2.3%.  For  the  TE21  mode,  the  efficiency 
drops  from  17 ~ 6. 85%  to  r;~3.76%  as  A y2/y0  increases 
to  2.3%.  The  efficiency  of  the  TM2!  mode  decreases  from 
r?  ~  2. 45  %  to  17  ~  1 . 27%  as  Ay2  /y 0  increases  to  1  % .  This 
is  more  rapid  than  for  the  TEm  and  TE21  modes,  and 
occurs  because  the  transition  to  the  thermal-beam  regime 
is  found  for  a  lower  value  of  the  momentum  spread  for  the 
TM21  mode. 

The  scaling  of  the  saturation  efficiency  with  beam 
current  for  the  TEoi,  TE21,  and  TM21  modes  is  shown  in 
Figs.  14-16,  respectively,  for  A y2  =0  and  Ay.  /y0  =  1%. 
On  the  basis  of  an  idealized  one-dimensional  model,29  it 
has  been  shown  that  the  saturation  efficiency  should  scale 
as  r)~iy}  for  frequencies  corresponding  to  peak  growth 
rates.  As  shown  in  Fig.  14,  this  type  of  scaling  law  is  ob¬ 
tained  for  the  TE0|  mode  over  a  range  of  currents  extend¬ 
ing  from  approximately  300-1000  A.  Observe  that  the 
dielectric  effect  of  the  electron  beam  on  the  waveguide 
mode  is  included  in  the  formulation,  and  the  bandwidth  of 
the  interaction  shifts  with  the  beam  current.  As  a  conse¬ 
quence,  the  increasing  divergence  between  the  simulation 
results  and  the  scaling  law  for  currents  below  300  A  is  at¬ 
tributed  to  a  shift  in  the  frequency  of  peak  growth  away 
from  co/ckw  =  11.3.  It  should  also  be  noted  that  the  curve 
for  A y./y0=  1%  increases  with  current  only  slightly  fas¬ 
ter  than  that  obtained  for  Ay2=0  and,  other  parameters 
being  equal,  the  1 1  '  scaling  law  seems  to  be  relatively  in- 
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FIG.  13.  Plot  of  the  saturation  efficiencies  of  the  TEm,  TE2i, 
and  TM;i  modes  vs  axial  energy  spread  for  ai/ckw  =  11.3. 


FIG.  14.  Graph  of  the  saturation  efficiency  of  the  TEoi  mode 
vs  beam  current  for  ai/ckw  =  11.3. 
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dependent  of  the  momentum  spread.  However,  one  point 
of  caution  ts  worth  noting  in  regard  to  these  conclusions. 
Specifically,  if  the  variation  in  current  were  accompanied 
by  variations  in  the  momentum  spread  and/or  beam  ra¬ 
dius  fas  might  be  expected  if  the  beam  were  held  a*,  con¬ 
stant  brightness),  then  the  scaling  law  might  differ.  The 
U/i  scaling  law  is  also  approximately  found  for  the  TF2: 


TMj!  Mtjdg  !s-8.8  om;  b-2.9  cm;  i/;4Kw«11.3;  P^-60  kW! 

E - -  - - 


.  t 


/ 


f  Vt  -  J  ■  M.V 

*  B0  -  1  6  am 

*  »„  •  m  »o 

»  »,l  sm  1 

.  N„  - J  j 

iw  m  m  3b'  Mo  mo  too  «o  mo 

l„.AI 

FIG.  ifi  Graph  of  the  saturation  efficiency  of  the  TM;i  mode 
vs  beam  current  for  w/cfc*  =  1 1.3. 


mode  for  currents  ranging  from  about  200-1000  A.  The 
slightly  larger  range  of  currents  over  which  the  scaling 
law  seems  to  hold  fox  the  TE2;  mode  is  attributed  10  ilie 
previously  mentioned  result  that  a/ckw^\\.l  is  closer 
to  the  frequency  of  peak  growth  late  for  this  mode  than 
for  the  TEoi  mode.  As  a  consequence,  this  frequency 
remains  in  the  vicinity  of  the  maximum  growth  rate  over 
a  broader  range  of  currents  for  the  TEji  mode.  Al¬ 
though  the  expected  scaling  law  is  approximately  ob¬ 
tained  for  the  TE01  and  TE;,  modes,  three-dimensional 
effects  involved  in  the  beam  propagation  and  the  cou¬ 
pling  between  the  beam  and  the  transverse  mode  struc¬ 
ture  do  modify  the  scaling.  This  is  most  evident  for  the 
TM21  mode  shown  in  Fig.  16  for  which  we  find  the 
efficiency  increasing  somewhat  faster  than  I,}  ■ 

We  now  address  the  issue  of  efficiency  enhancement  by 
means  of  a  tapered  wiggler  field,  and  concentrate  on  the 
TEo;  mode.  The  fundamental  theory  of  the  efficiency 
enhancement  mechanism  has  been  amply  discussed  in  the 
literature. 51 -3f  In  order  to  understand  the  physical  basis 
of  the  process,  we  observe  that  the  wave-particle  reso¬ 
nance  condition  which  gives  nse  to  amplification  is  ex¬ 
tremely  sensitive  to  the  axial  electron  velocity.  However, 
the  axial  velocity  decreases  as  the  wave  is  amplified  until, 
ultimately,  the  resonance  condition  is  broken  and  the 
amplification  ceases.  The  purpose  of  the  tapered  wiggler 
field  is  to  reduce  the  transverse  velocity  induced  by  the 
wiggler  which,  in  turn,  results  in  an  axial  acceleration  that 
maintains  the  resonance  condition.  Thus,  the  tapered 
field  is  a  means  of  “tapping”  the  transverse  kinetic  energy 
of  the  Deam.  This  can  be  accomplished  by  a  tapering  of 
either  the  wiggler  amplitude  or  period;  however,  in  the 
present  work  we  shall  confine  ourselves  to  a  tapered  am¬ 
plitude  configuration.  The  results  of  the  simulation  are  in 
qualitative  agreement  with  those  obtained  for  a  helical 
wiggler  configuration,24  and  indicate  that  the  efficiency 
enhancement  is  extremely  sensitive  both  to  the  point  at 
which  the  taper  is  begun  and  to  the  slope  of  the  taper. 
The  optimal  position  at  which  to  begin  the  taper  is,  typi¬ 
cally,  at  a  point  shortly  prior  to  saturation  (for  the  unta¬ 
pered  system)  which  corresponds  to  the  trapping  of  the 
bulk  of  the  electron  beam  in  the  pondcromotive  potential 
formed  by  the  beating  of  the  wiggler  and  radiation  fields. 

For  the  case  corresponding  to  the  TE01  mode  shown  in 
Fig.  2,  the  optimal  point  at  which  to  begin  the  taper 
Occurs  for  —  1 10.  In  order  to  accelerate  the  beam  in 
the  axial  direction  the  wiggler  field  must  be  decreased, 
and  the  evolution  of  the  wave  power  versus  axial  position 
is  shown  in  Fig  1?  for  a  taper  of  *  —0.007.  The 
efficiency  enhancement  (as  measured  by  the  growth  in  the 
wave  power!  continues  for  as  long  as  the  taper  is  main¬ 
tained,  and  a  maximum  efficiency  of  t7mM=I34.0% 's  ob¬ 
tained  if  the  wiggler  field  is  tapered  to  zero  at  253 
(i.e.,  3.95  mi.  The  effect  of  different  choices  for  the 
wiggler  taper  is  shown  in  Fig.  18,  in  which  we  plot  the 
maximum  efficiency  found  by  tapering  the  wiggler  field  to 
zero  versus  c„ .  As  shown  in  the  figure,  optimal  results 
are  found  for  e w  2=  —0.002  at  which  point  and 

the  interaction  region  extends  to  <cwz=r610  (i.e„  9.5  m). 
The  variation  in  the  relative  phase  for  a  tapered  wiggler 
interaction  is  shown  in  Fig.  19,  in  which  we  plot  the  reia- 
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FIG.  1 7.  Plot  of  the  evolution  of  the  TE0,  mode  for  a  tapered 
wiggler  interaction  with  e„,  =  —0.007  and  k.rz o=  1 10. 

tive  phase  versus  axial  position  for  eu  = —0.007,  and 
a)/ck,x  =  10.7,  11.0,  and  11.3.  Note  that  the  start  taper 
point  is  chosen  to  be  the  optimum  value  for  each  frequen¬ 
cy.  and  is  indicated  in  the  figure  by  the  airow.  Three 
features  are  readily  apparent  from  the  figure.  The  first  is 
that  the  variation  in  the  relative  phase  subsequent  to  the 
start-taper  point  is  approximately  the  same  for  each  of  the 
three  frequencies  shown,  and  the  spacing  between  the 
curves  remains  approximately  constant.  Second,  while 
the  phase  variation  within  the  tapered  wiggler  region  ini¬ 
tially  increases  immediately  after  the  start-taper  position, 
the  phase  variation  appears  to  saturate  and  remain  rela¬ 
tively  constant  over  an  extended  interaction  length. 
Third,  the  oscillations  at  one-half  the  wiggler  appear  to 
decrease  in  amplitude  over  the  course  of  the  tapered 
wiggler  region. 
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FIG.  18.  Graph  of  the  maximum  efficiency  of  the  TE„i  mode 
vs  wiggler  taper  (f». 1. 
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FIG.  19.  Plot  of  the  evolution  of  the  relative  phase  during  the 
course  of  the  tapered  wiggler  interaction  for  io/ck^.  =  10.7,  11.0, 
and  1 1.3. 


The  effect  of  an  initial  beam  momentum  spread  on  the 
tapered  wiggler-efficiency  enhancement  process  has  also 
been  investigated,  and  the  results  are  summarized  in  Fig. 
20  in  which  we  plot  the  maximum  realizable  efficiency  (if 
the  wiggler  field  is  tapered  to  zero)  versus  Ay.  /y0.  The 
maximum  efficiency  at  a  fixed  start-taper  point 
(fc.,z0  =  110)  chosen  to  correspond  to  the  optimum  posi¬ 
tion  for  Ay.=0  is  shown  in  Fig.  20(a)  versus  the  axial 
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FIG.  20.  Illustration  of  the  effect  of  beam  momenlum 
spread  on  the  tapered  wiggler  interaction  at  (a)  fixed  z0.  and  ' bi 
variable  z0-  Observe  that  in  (b)  the  optimal  start-taper  points 
correspond  to  <:,lz0  =  llO  for  Ay,/yo<l%,  ku.20  =  ll7  for 
Ay,/)'o=  15%,  and  ku.z0  =  124  for  Ay,  /y0  =  2%. 
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momentum  spread.  As  is  evident  from  the  figure,  the 
efficiency-enhancement  process  is  unaffected  by  the 
momentum  spread  for  d'/J/yc,C>  1%,  but  decreases  rapid¬ 
ly  for  axial  energy  spreads  above  this  value.  There  are 
two  principal  reasons  for  this  decrease  in  The  first 

is  that  the  phase-trapping  mechanism  becomes  less 
effective  as  the  momentum  spread  increases  because  a  pro¬ 
portionally  greater  fraction  of  the  beam  remains  outside 
the  trapped  orbit  region  of  the  axial  phase  space.  The 
second  reason  is  that  the  increase  in  the  momentum 
spread  results  in  a  decrease  in  the  growth  rate  ana  a 
longer  distance  to  saturation.  Because  of  this,  the  op¬ 
timum  stan-taper  point  is  an  increasing  function  of  the 
momentum  spread.  Thus,  if  we  determine  the  optimum 
start-taper  point  versus  Ay,,  then  the  maximum  efficiency 
can  be  expected  to  decrease  less  rapidly  with  axial  energy- 
spread  than  is  shown  in  Fig.  20(a).  This  is  indeed  the 
case,  as  shown  in  Fig.  20(b)  in  which  we  plot  t?,n„  versus 
Ay,  for  the  optimum  stan-taper  points.  Note  that  for 
Ay,/y0Sl%  this  point  coincides  with  that  for  Ay, -0 
Tor  the  current  choice  of  parameters).  As  seen  in  the 
figure,  7jfrM~26.6‘%)  for  Ay,/y3  =  2%  and  a  start-taper 
point  of  ku2c  =  124.  This  contrasts  with  an  rjinu  1 1 . 8 9c 
'and  Ay,  /y0=2%i  when  Ic,„z0  =  110.  As  a  consequence, 
we  conclude  that  although  the  untapered  efficiency  de¬ 
creases  relatively  quickly  with  increasing  momentum 
spread,  the  tapered-wiggler  interaction  can  accept  a  small 
level  of  momentum  spread  without  suffering  a  significant 
degradation  in  the  interaction  efficiency.  However,  the 
upper  limit  on  the  allowable  momentum  spread  must  be 
determined  by  the  particular  choice  of  experimental  pa¬ 
rameters,  and  the  value  of  Ay,/y<)Sl%  should  net  be 
construed  to  be  a  general  result. 

IV.  SUMMARY  AND  DISCUSSION 

In  this  paper  a  fully  self-consistent  nonlinear  theory  and 
numerical  simulation  has  been  developed  for  the  FEL 
amplifier  m  three  dimensions.  The  particular  config¬ 
uration  of  interest  consists  of  a  cylindricaliy  symmetric 
electron  beam  of  arbitrary  cross  section  ion  entry  at  2  =0) 
injected  into  a  loss-free  rectangular  waveguide  in  the  pres¬ 
ence  of  a  hnearly  polarized  wiggler  magnetic  field.  The 
wiggler-fleld  model  is  that  generated  by  a  magnet  stack  in 
which  the  individual  magnets  have  parabolically  tapered 
pole  pieces  in  order  to  provide  for  electron  focusing  in  the 
plane  of  the  bulk  wiggler-induced  motion.  In  addition, 
the  adiabatic  injection  of  the  electron  beam  is  modeled  by 
the  inclusion  of  an  initial  taper  of  the  wiggler  amplitude. 
The  system  of  equations  derived  is  a  generalization  of  a 
previously  described  analysis  of  a  helical  wiggler 
cylindrical  waveguide  configuration”  ~v  which  treats 
the  self-consistent  evolution  of  the  trajectories  of  an  en¬ 
semble  of  electrons  and  the  radiation  fields.  The  analysis 
includes  both  the  T£  and  TM  modes,  and  includes  the 
overlap  of  the  transverse  mode  structure  and  the  electron 
beam  (i.e„  the  filling  factor)  in  a  self-consistent  way. 
Space-charge  fields  have  been  neglected;  hence,  the 
analysis  is  restricted  to  the  mgh-gain  Compton  (strong- 
pump)  regime  of  operation.  Since  the  problem  of  in¬ 
terest  is  the  FEL  amplifier,  only  single-frequency  propa¬ 


gation  is  considered,  which  permits  an  average  over  a 
wave  period  to  be  performed  that  eliminates  the  fast- 
time-scale  phenomena  formulation.  This  results  in  a 
great  increase  in  computational  efficiency  over  a  full- 
scale  particle-in-all  simulation  code,  and  allows  the  ap 
plication  of  the  technique  to  short  ti.e.,  optical)  wave¬ 
lengths  given  an  appropriate  mode  structure. 

The  electron  trajectories  are  integrated  using  the  com¬ 
plete  Lorentz-force  equations,  so  we  are  able  to  study  the 
detailed  orbital  dynamics  in  the  combined  wiggler- 
radiation  field  structure.  The  overall  bulk  motion  of  the 
electron  beam  exhibits  a  dominant  oscillation  at  the 
wiggler  period,  as  well  as  a  slow-time-scalc  betatron  oscil¬ 
lation  due  to  the  transverse  inhomogeneity  in  the  wiggler 
field.  The  dynamics  of  the  adiabatic  injection  of  the  elec¬ 
tron  beam  were  studied  by  means  of  a  comparison  of  the 
saturation  efficiency  with  the  length  of  the  entry  taper  re¬ 
gion.  The  results  indicate  that  the  saturation  efficiency  in¬ 
creases  relatively  quickly  with  the  length  of  the  entry- 
taper  region  for  \u.  6,  and  more  slowly  thereafter.  The 

reason  for  this  is  that  the  injection  process  itself  in  a  real¬ 
istic  (i.e.,  three-dimensional)  wiggler  field  introduces  a 
effective  momentum  spread  on  the  beam  which  decreases 
as  the  axial  wiggler  gradient  becomes  more  gradual.  For 
practical  purposes  .V,„  £  S  appears  to  be  an  adequate 
compromise  between  the  minimization  of  the  overall 
wiggler  length  and  the  maximization  of  the  interaction 
efficiency,  and  this  is  the  regime  in  which  many  mi¬ 
crowave  FEL’s  have  been  operated/-1-10"”  In  addition, 
the  effect  of  a  tapered  wiggler  amplitude  on  the  enhance¬ 
ment  Of  the  interaction  efficiency  can  be  included  in  a 
straightforward  manner. 

The  numerical  analysis  has  been  performed  for  the  case 
of  a  35-GHz  FEL  amplifier  which  employs  a  3.3-MeV, 
800-A  electron  beam  with  an  initial  radius  of  1.0  cm. 
The  beam  propagates  through  a  rectangular  waveguide 
with  dimensions  a  =9.8  cm  and  6=2.9  cm  in  the  pres¬ 
ence  of  a  wiggler  field  with  a  3.72-kG  amplitude  and  9.8- 
cm  period.  Three  distinct  waveguide  modes  are  found  to 
bs  amplified:  ihe  TE0|,  TE21,  and  TMj|  modes.  For  this 
choice  of  frequency,  the  TE0;  mode  exhibits  the  highest 
(untapered;  efficiency,  while  the  TE;i  mode  has  the 
highest  growth  rate.  The  wave-panicle  coupling  for  the 
TMj;  mode  is  the  weakest  of  the  three  modes  and  has  the 
lowest  growth  rate  and  efficiency  as  well  as  the  narrowest 
bandwidth.  The  effect  of  an  initial  momentum  spread  is 
investigated  for  axial  energy  spreads  Ay,,  y0S2%  which, 
for  this  choice  of  parameters,  is  well  within  the  cold-beam 
regime.  Results  indicate  that  over  this  range  of  Ay,  the 
efficiency  decreases  in  an  approximately  linear  fashion 
with  increasing  axial  energy  spread.  This  is  in  substantial 
agreement  with  the  results  obtained  for  a  helical  wiggler 
configuration;”  however,  it  should  be  remarked  that  the 
same  distribution  has  been  used  for  both  the  planar  and 
helical  wiggler  configurations.  Thus,  while  the  planar  and 
helical  wigglers  behave  in  substantially  the  same  manner, 
the  detailed  scaling  of  the  efficiency  with  the  momentum 
spread  can  be  expected  to  vary  with  the  detailed  choice  of 
distribution. 

The  phase  variation  of  the  planar  wiggler  configuration 
discussed  herein  is  also  in  qualitative  agreement  with  that 
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found  for  helical  wiggler  configurations.26  The  principal 
difference  is  that  an  oscillation  at  one-half  the  wiggler 
period  is  superimposed  on  the  bulk  variation  in  the  phase 
due  to  details  of  the  wave-particle  coupling  with  a  planar 
wiggler.  Specifically,  for  the  gain  band  associated  with 
the  upper  (high-frequency)  intersection  between  the 
beam-resonance  line  and  the  waveguide-dispersion  curve, 
we  observe  that  the  bulk  variation  of  the  relative  phase  de¬ 
creases  with  axial  position  up  to  a  point  short  of  that  at 
which  the  power  saturates  for  frequencies  at  the  low  end 
of  the  band.  As  the  frequency  increases,  the  phase  varia¬ 
tion  decreases  until  a  critical  frequency  is  reached  for 
which  the  phase  remains  relatively  constant  over  the 
course  of  the  interaction.  This  critical  frequency  has  been 
found  to  occur  at  approximately  10%  below  the  frequen¬ 
cy  of  peak  growth  rate  for  all  parametric  cases  studied  for 
both  the  helical  and  planar  wigglers.  Above  this  critical 
frequency  the  average  relative  phase  tends  to  increase  with 
axial  position.  In  view  of  the  high  power  potential  of  the 
FEL  amplifier,  applications  such  as  microwave  sources  for 
the  next  generation  of  radio-frequency  electron  accelera¬ 
tors  are  natural  considerations  and  the  question  of  the 
phase  stability  of  these  devices  against  fluctuations  in  the 
electron-beam  voltage  is  of  importance.  Again,  we  find 
qualitative  agreement  on  this  issue  between  simulations  of 
helical  and  planar  wiggler  configurations.  The  results  in¬ 
dicate  a  much  poorer  phase  stability  than  the  current  gen¬ 
eration  of  SLAC  klystrons,  and  we  conclude  that  applica¬ 
tions  of  FEL  amplifiers  which  require  an  extremely 
phase-stable  microwave  source  will  also  require  an  elec¬ 
tron  beam  with  a  very  low  level  of  voltage  fluctuations. 

The  enhancement  of  the  interaction  efficiency  by  means 
of  a  tapered  wiggler  amplitude  shows  maximum 
efficiencies  of  the  order  of  35-45  %  are  possible  for  this 
choice  of  parameters.  This  brings  maximum  power  levels 
into  the  GW  range.  The  simulation  results  also  indicate 
that  the  tapered  wiggler-efficiency-enhancement  mecha¬ 
nism  is  relatively  less  sensitive  to  the  effect  of  momentum 
spread  than  the  uniform-wiggler  case,  and  no  degradation 
in  the  maximum  efficiency  is  found  for  ^y,/yQS  1%  in 
the  present  case.  Examination  of  the  phase  variation  dur¬ 
ing  the  tapered  wiggler  interaction  shows  similar  results 
over  the  entire  gain  band  (Fig.  19),  and  there  is  no  de¬ 
crease  in  the  phase  separation  for  the  various  frequencies 
observed  in  the  tapered  wiggler  region.  Since  variations  in 
the  frequency  at  fixed  beam  voltage  are  fundamentally 
equivalent  to  a  variation  in  the  voltage  at  fixed  frequency 
(i.e.,  both  processes  sweep  through  the  gain  band),  we  find 
that  our  conclusion  regarding  the  phase  stability  of  FEL 
amplifiers  is  unaltered  for  the  case  of  tapered  wiggler  in¬ 
teractions. 

Indeed,  the  bulk  features  of  the  FEL  interaction  mecha¬ 
nism  are  found  to  be  in  substantial  agreement  for  both 
helical  and  planar  wiggler  configurations.  The  most  evi¬ 
dent  distinction  is  the  rapid  oscillation  in  both  the  power 
and  relative  phase  which  appears  for  planar  wigglers  at 
half  the  wiggler  period.  While  the  bulk  wave 
amplification  is  unaffected  by  this  oscillation,  it  may  intro¬ 
duce  an  uncertainty  in  the  measurement  of  the  gain 
and/or  power  from  planar  wiggler  configurations  which  is 
of  the  order  of  10%  (apart  from  systematic  instrumental 


errors).  In  addition,  the  relative  phase  also  exhibits  this 
oscillation,  which  indicates  a  periodic  modulation  in  wave 
refraction  also  occurs.  It  is,  therefore,  an  open  question 
whether  this  may  affect  the  focusing  (i.e.,  optical  guiding) 
of  the  radiation,  and  this  will  be  addressed  in  a  future 
work  by  the  inclusion  of  multiple  modes  in  the  formula¬ 
tion.  However,  this  question  may  be  moot  for  tapered 
wiggler  configurations  since  the  oscillation  appears  to  be 
attenuated  (Fig.  19). 

Although  the  configuration  described  in  this  paper  does 
not  precisely  coincide  with  the  experiment  conducted  by 
Orzechowski  and  co-workers,913  the  parameters  chosen 
for  the  numerical  analysis  coincide  with  those  of  the  ex¬ 
periment  and  it  is  useful  to  compare  the  simulation  with 
the  experiment.  The  fundamental  differences  between  the 
experimental  configuration  and  the  analytical  model  are 
that  in  the  experiment  (1)  the  beam  was  injected  into  the 
wiggler  by  means  of  an  entry  taper  region  one  wiggler 
period  in  length,  (2)  a  quadrupole  field  was  used  to  pro¬ 
vide  electron  focusing  instead  of  parabolically  shaped  pole 
pieces,  and  (3)  in  the  tapered  wiggler  experiment13  the  am¬ 
plitude  was  tapered  nonlinearly.  A  50-kW  34.6-GHz  (i.e., 
a>/ckw  =  1 1.3)  magnetron  was  used  to  drive  the  FEL 
amplifier,  and  the  signal  was  injected  in  such  a  way  as  to 
couple  primarily  to  the  TE01  mode.  As  a  result,  the  TE01 
mode  was  the  predominant  component  of  the  output  sig¬ 
nal;  however,  significant  power  levels  were  also  detected 
in  the  TE;i  and  TM21  modes.  The  detailed  experimental 
parameters  we  choose  to  compare  with  the  experiment  in¬ 
volve  a  3.5-MeV,  850-A  electron  beam  used  in  the  tapered 
wiggler  experiment.  Although  the  beam  is  thought  to  ex¬ 
hibit  a  slightly  elliptical  cross  section,  the  average  of  the 
semimajor  and  semiminor  radii  is  about  1.0  cm  and  we 
choose  this  value  for  the  beam  radius  in  the  numerical 
simulation.  The  waveguide  and  wiggler  parameters  are 
those  used  previously  in  Sec.  Ill,  and  we  note  that  since 
the  fringing  fields  are  not  included  in  our  model  of  the  en¬ 
try  taper  region  the  choice  of  .Vu.  =5  is  a  compromise. 
The  evolution  of  the  wave  power  in  the  TE0t  mode  as  a 
function  of  axial  position  for  these  parameters  is  shown  in 
Fig.  21  for  a  uniform  wiggler  (6^=0),  and  for 
e„.  =  —  0.007.  The  results  are  similar  to  those  described 
in  the  preceding  section  for  an  800-A  beam,  and  the  best 
fit  with  the  experiment  is  found  for  an  axial  energy  spread 
of  Ayz/y0^  As  shown  in  the  figure,  the  untapered 
wiggler  results  give  a  saturated  power  of  194  MW  over  a 
uniform  wiggler-interaction  region  (i.e.,  Zu,  —  )  of  1.4 

m  in  length.  This  compares  well  with  the  experimental 
measurement  of  180  MW  over  an  interaction  length  of  1.3 
m,  and  the  discrepancies  are  within  the  experimental  un¬ 
certainty.  We  remark,  however,  that  the  simulation  re¬ 
sults  may  be  affected  by  as  much  as  20%  by  the  choice  of 
the  length  of  the  entry  taper  region.  It  is  more  difficult  to 
compare  the  simulation  (with  a  linear  wiggler  taper)  and 
the  experiment  (with  a  nonlinear  taper)  in  such  a  detailed 
fashion.  It  is  certainly  reasonable  to  expect  that  a  non¬ 
linear  taper  might  be  more  efficient  than  a  linear  one  in 
the  sense  that  the  efficiency  enhancement  can  be  accom¬ 
plished  over  a  shorter  interaction  region.  Thus,  in  the 
comparison  of  the  tapered  wiggler  results  we  interpret  the 
simulation  as  an  indication  of  limits  on  the  energy  extrac- 
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FIG.  21  Plot  showing  the  evolution  of  the  TE0i  mode  vs  axi¬ 
al  position  for  a  uniform  and  tapered  higgler  interaction. 


non.  With  this  in  mind  we  observe  that  the  optimal 
start-taper  point  is  at  k,^z0  &  1 1 1  (slightly  more  than  one 
wiggler  period  short  of  the  saturation  point  for  a  uniform 
wiggleri.  We  choose  the  slope  of  the  taper  to  correspond 
roughly  with  the  average  taper  in  the  experiment,  and  find 
a  maximum  efficiency  of  about  35. 6%  at  a  power 

level  of  1.06  GW'.  This  differs  by  only  about  %%  from 
that  found  in  the  experiment,  and  we  interpret  this  as 
good  agreement  within  the  uncertainty  introduced  by  the 
differences  in  the  slope  of  the  taper  We  conclude,  there¬ 
fore.  that  the  simulation  is  consistent  with  the  interpreta¬ 
tion  of  a  1%  axial  beam  energy  spread,  which  is  support¬ 
ed  by  an  electron-spectrometer  measurement  indicating  an 
upper  limit  on  the  energy  spread  of  2%.’°  Furthermore, 
■n  view  of  the  relative  insensitivity  of  the  tapered  wiggler 
interaction  (see  Fig.  20:.  we  conclude  that  the  perfor¬ 
mance  of  the  tapered  wiggler  experiment  would  not  be 
markedly  improved  by  a  further  improvement  in  beam 
quality  (i.e.,  a  decrease  in  the  momentum  spread). 

It  is  useful  to  compare  the  analysis  and  simulation  de¬ 
scribed  herein  with  other  models  in  order  to  assess  the  rel¬ 
ative  merits  of  the  various  formulations  of  the  problem 
To  this  end,  we  consider  the  simulation  code  FRED  in  use 
at  the  Lawrence  Livermore  National  Laboratory28,40  be¬ 
cause  this  code  has  been  used  extensively  to  model  and  in¬ 


terpret  the  experiment.  It  is  important  to  recognize  that 
simulations!  modeh  and  computer  codes  undergo  lapid 
development  which  renders  an  accurate  description 
difficult;  hence,  1  will  confine  my  remarks  to  fundamental 
properties  which  include  the  particle  dynamics  and  ihe 
treatment  of  the  radiation  field.  The  particle  dynamics  in 
FRED  are  treated  m  a  reduced  form  by  averaging  the  orbit 
equations  over  a  wiggler  period.  The  transverse  motion  is 
described  in  terms  of  the  bulk  wiggler  and  betatron  oscil¬ 
lations  by  means  of  an  analytical  approximation,  and  the 
self-consistent  effect  of  the  radiation  field  on  the  trans¬ 
verse  motion  is  not  included.  The  dynamics  in  the  axial 
direction  are  handled  by  integration  of  equations  for  the 
particle  phase  and  energy.  This  is  the  most  common  ap¬ 
proach  used  in  the  simulation  of  free-electron  lasers  and, 
as  shown  by  the  comparison  between  FRED  and  the  exper¬ 
imental  results,’’ :i  it  works  well.  However,  certain  effects 
are  excluded  from  this  formuiation.  First,  we  observe 
that  the  rapid  oscillation  at  the  second  wiggler  harmonic 
is  lost  entirely  due  to  the  averaging  of  the  orbit  equations 
over  a  wiggler  period.  Second,  the  injection  of  the  elec¬ 
tron  beam  into  the  wiggler  is  excluded  from  the  simula¬ 
tion,  and  must  be  described  explicitly.  Third,  because  the 
self-consistent  effect  of  the  radiation  field  on  the  trans¬ 
verse  electron  motion  is  net  included,  the  initial  transients 
associated  with  the  injection  of  the  radiation  into  the 
amplifier  are  not  properly  handled.  As  a  result,  launching 
losses  cannot  be  described  by  FRED,  and  the  wave  power 
must  be  initialized  at  a  lower  value  than  that  actually  in¬ 
jected  into  the  amplifier.0  Thus,  it  is  our  opinion  that  the 
approach  to  the  particle  dynamics  described  in  this  paper 
is  superior  to  that  employed  in  FRED.  The  advantage  of 
FRED  over  the  single-mode  analysis  we  have  described  is 
that  the  radiation  field  is  treated  by  means  of  a  field  solver 
which  implicitly  handles  multiple  modes  and,  thus,  de¬ 
scribes  the  optical  guiding  (i.e.,  focusing)  of  the  radiation 
due  to  the  electron-beam  interaction,  However,  optical 
guiding  can  be  treated  via  a  multirnoae  formulation  of  the 
present  analysis  as  well.  Indeed,  such  an  analysis  is  now 
under  study,  and  the  results  will  be  presented  in  the  forth¬ 
coming  work. 
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A  nonlinear  formulation  of  the  free-electron-laser  amplifier  with  a  linearly  polarized  wiggler 
magnetic  field  is  used  to  study  harmonic  generation.  Substantial  emission  is  found  to  occur  at  the 
harmonics  for  a  cold  beam;  however,  the  harmonics  are  far  more  sensitive  to  beam  thermal  effects 
than  is  the  fundamental. 


The  question  of  harmonic  generation  in  free-electron 
lasers  (FEL’s)  has  become  increasingly  important  in  re¬ 
cent  years.  The  operation  of  an  FEL  oscillator  at  visible 
wavelengths  driven  by  the  ACO  (les  Anneaux  de  Col¬ 
lisions  de  l’Accelerateur  Lineaire  d’Orsay)  storage  ring1 
has  shown  that  the  generation  of  harmonics  at  uv  wave¬ 
lengths  can  limit  performance  due  to  severe  degradation 
of  the  optical  system.  In  contrast,  there  is  interest  in  the 
generation  of  harmonics  in  an  effort  to  achieve  ever 
shorter  wavelengths  without  the  necessity  of  ever  higher 
electron  beam  energies.  In  order  to  address  these  ques¬ 
tions,  we  have  developed  a  fully  three-dimensional  for¬ 
mulation  of  the  FEL  amplifier  for  a  configuration  which 
consists  of  an  electron  bsam  propagating  through  a  rec¬ 
tangular  waveguide  in  the  presence  of  a  linearly  polar¬ 
ized  wiggler  magnetic  field.  A  set  of  nonlinear 
differential  equations  is  derived  which  self-consistently 
describes  the  evolution  of  both  an  ensemble  of  electrons 
and  the  electromagnetic  fields.  The  issue  of  harmonic 
generation  [i.e.,  co^ik  +lkw  )v:  for  l>  1]  has  been  ad¬ 
dressed  previously  for  the  case  of  an  electron  beam  prop¬ 
agating  through  a  cylindrical  waveguide  in  the  presence 
_ l 


of  a  helical  wiggler  field,2  and  emission  was  shown  to 
occur  at  all  harmonics  of  the  resonant  frequency.  How¬ 
ever,  due  to  the  modulation  of  the  axial  electron  velocity 
in  linearly  polarized  wiggler  fields,  only  the  odd  harmon¬ 
ics  are  excited.  In  this  paper  we  study  both  the  linear 
growth  rates  and  nonlinear  saturation  efficiencies  at  the 
odd  harmonics,  as  well  as  the  effect  of  thermal  spread  on 
harmonic  generation. 

The  configuration  we  consider  is  that  of  a  relativistic 
electron  beam  propagating  through  a  loss-free  rectangu¬ 
lar  waveguide  in  the  presence  of  a  linearly  polarized 
wiggler  magnetic  field.  The  wiggler-field  model  we  em¬ 
ploy  is  that  encountered  when  the  individual  magnets  in 
the  wiggler  have  parabolically  tapered  pole  faces.  The 
wiggler  field  generated  in  this  way  provides  for  electron 
beam  focusing  in  the  plane  of  the  principal  wiggler 
motion,  and  was  first  employed  experimentally  by  Phil¬ 
lips.3  A  detailed  analysis  of  the  magnetic  field  produced 
by  a  wiggler  with  parabolically  tapered  pole  pieces  was 
undertaken  by  Scharlemann,4  who  showed  that  the 
wiggler  field  is  of  the  form 
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where  Bw  denotes  the  wiggler  amplitude  and 
kw  (  =  2ir/kw)  is  the  wiggler  wave  number.  In  addition, 
we  model  the  injection  of  an  electron  beam  into  the 
wiggler  by  allowing  the  wiggler  amplitude  to  increase 
adiabatically  from  zero  to  a  constant  level  over  Nw 
wiggler  periods  as  follows: 


Bjz)  = 


Bw  sin2 


kwz 


4  N, 

Bw,  z  >Swkw  . 


,  0<z^,\wkw 


(2) 


Note,  however,  that  fringing  fields  associated  with  the 
tapered  wiggler  amplitude  are  neglected,  and  this  repre¬ 
sentation  (for  z  <Nwkw)  is  valid  only  for  large  Sw.  For 
the  purposes  of  the  present  analysts,  we  choose  Nw  =  10. 

The  boundary  conditions  at  the  waveguide  wall  may 
be  satisfied  by  expanding  the  vector  potential  in  terms  of 
the  orthogonal  basis  functions  of  the  TE  modes,  and  we 
write  the  vector  potential  of  the  radiation  field  in  the 
form 

6  A(i,f)  =  2'  6/4MB(z)ems(x,y)cosa  ,  (3) 

m  «0 

*  .0 
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where  for  frequency  a>  and  wave  number  k  ( z ) 

a=  f:dz'k(z‘)-ot  , 

J  o 

J'  indicates  that  m  and  n  are  not  both  zero,  and 
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is  the  polarization  vector.  In  this  representation  the 
waveguide  is  assumed  to  be  centered  at  the  origin  and 
bounded  by  —a/2<^x<a/2  and  —  b/2<y  <b  /2.  As  a 
consequence,  X=x+a/ 2,  Y  =y  +  b/ 2,  and 


km„  =ir 


1/2 


(6) 


denotes  the  cutoff  wave  number.  It  is  implicitly  assumed 
that  both  the  mode  amplitudes  8Am„(z)  and  the  wave 
number  k  (z)  vary  slowly  over  a  wave  period. 

A  detailed  derivation  of  the  equations  which  describe 
the  evolution  of  the  mode  amplitudes  and  wave  number 
has  been  presented5  for  both  the  TE  and  TM  modes; 
however,  the  method  follows  that  employed  for  a  helical 
wiggler.2  Here  we  note  that  substitution  of  the  fields  (3) 
into  Maxwell’s  equations  yields 
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upon  orthogonalization  in  x  and  y,  and  an  average  over 
the  wave  period.  In  Eqs.  (7)  and  18), 
8amn=e8Amn/mc2,  ojb  is  the  ambient  plasma  frequen¬ 
cy,  0zO  =  vl0/c  (where  vz0  is  the  initial  axial  beam  veloci¬ 
ty),  v  is  the  instantaneous  velocity,  and  Fmn  =4  (8)  when 
either  m  =  0  or  n  =  0  (m# 0  and  n^O).  The  averaging 
operator  ((  •  ■  •  )}  describes  the  source  current  and  in¬ 
cludes  the  effect  of  an  initial  (i.e.,  at  z=0)  thermal 
spread  in  the  beam  by  means  of  the  distribution  function 

F0(p0)=A  exp[  —  (pl0  —p0  )2 / A/?*2 ] 

x8(po-p2LO-P?oWP:o)  .  l9) 


where  p 0  and  A pt  describe  the  initial  bulk  momentum 
and  the  momentum  spread  of  the  beam,  0(x)  is  the 
Heaviside  function,  and  the  normalization  constant  is 

A  ~  l71-  f*°  dPzoexp[-{Pio-Po)2/hp:]  j  .  (10) 

This  distribution  describes  an  initially  monoenergetic 
beam  with  a  pitch-angle  spread.  The  average  is  over  the 
initial  conditions  and  defined  as 


<  ••  )  =  J^ddo  dp.0-^- sxp[-(p.0-p0)2/±p:] 

lirab  Jn  n « 


P  o 


JT  daocr^cto)  fa2dx0  'dy0crLx0,y0H 


where  a0  (  =  -<ur0  and  f0  is  the  time  the  electron  crosses 
the  z  =0  plane)  is  the  initial  phase,  i£0  =  tan  "  '(pyo/pio), 
and  ct  (a0 1  and  c7L(x0,y0)  describe  the  electron  distribu¬ 
tions  in  initial  phase  and  initial  cross  section.  For  sim¬ 
plicity,  we  assume  the  beam  to  be  uniformly  distributed 
in  phase  ( <r l(  =  1 )  and  cross  section  for  an  initially  cylin¬ 
drical  beam  (i.e.,  <72=1  for  r0<Rb).  In  order  to  com¬ 
plete  the  formulation,  the  electron  trajectories  are  fol¬ 
lowed  by  means  of  the  Lorentz  force  equations  in  the 
combined  wiggler  and  radiation  fields.  It  should  be 
remarked  that  although  the  fast-time-scale  behavior  of 
Maxwell’s  equations  have  been  eliminated  by  means  of 
the  average  over  a  wave  period,  no  average  of  the 
Lorentz  force  equations  is  undertaken.  This  is  impor¬ 
tant  in  order  to  resolve  the  interaction  at  higher  har¬ 
monics. 

Equations  (7)  and  (8)  can  be  reduced  to  a  set  of  three 
first-order  differential  equations  describing  the  evolution 
of  the  amplitude  8amnf  the  growth  rate  Tm, 
(  =d  \nbamn/dz),  and  the  wave  number.  Hence,  includ¬ 


ing  the  orbit  equations,  we  solve  an  initial-value  problem 
composed  of  a  set  of  6,V-t-3  coupled  nonlinear 
differential  equations,  where  .V  is  the  number  of  parti¬ 
cles.  The  field  quantities  are  initialized  such  that 
8am„iz  =0)  describes  an  arbitrary  input  power  level. 
T^,(z  =0)=0,  and  kiz  =0)  —  {co2/c2 -k1„„)wl.  Thus, 
while  the  initial  state  describes  a  vacuum  mode,  the  sys¬ 
tem  evolves  into  a  dielectrically  loaded  waveguide  mode 
in  the  presence  of  the  beam. 

The  particular  example  we  consider  involves  a  3.3- 
MeV,  100- A  electron  beam  with  an  initial  radius  of  0.2 
cm  propagating  through  a  rectangular  waveguide  (a  =  10 
cm,  6=3  cm)  in  the  presence  of  a  wiggler  with  a  4.2-kG 
amplitude  and  a  9.8-cm  period.  A  spectrum  of  the 
efficiency,  rj,  of  the  TEot  mode  is  shown  in  Fig.  1  as  a 
function  of  frequency  for  the  fundamental  through  the 
ninth  harmonic  in  the  absence  of  any  thermal  spread 
(i.e.,  Ap,=0).  The  doublet  at  the  fundamental  corre¬ 
sponds  to  the  upper  and  lower  intersections  between  the 
waveguide  dispersion  curve  and  the  beam-resonance  line 


3220 


H.  P.  FREUND,  C.  L.  CHANG,  AND  H.  BLUEM 


36 


TE01  Mode  (a  =  10  cm;  b  =  3  cm) 


FIG.  t.  Plot  of  the  efficiency  vs  frequency  at  the  fundamen¬ 
tal,  third,  fifth,  seventh,  and  ninth  harmonics. 


+kw  )u,  ].  It  is  evident  from  the  figure  that  there 
is  a  sharp  drop  in  the  efficiency  between  the  fundamental 
and  the  third  harmonic.  However,  the  subsequent  de¬ 
crease  of  the  efficiency  with  harmonic  number  is  rather 
slow,  and  substantial  power  may  be  found  at  the  higher 
harmonics.  In  addition,  the  decrease  in  the  growth  rates 
with  harmonic  number  is  relatively  slow.  B;cau«e  cf 
variations  in  the  launching  loss  associated  with  the  vari¬ 
ous  harmonics,  as  well  as  oscillations  in  the  growth  rate 
over  the  course  of  the  interaction,  it  is  difficult  to  ob'ain 
estimates  of  the  average  growth  rates.  However,  the  dis¬ 
tance  to  saturation  lzsal )  for  a  fixed  input  signal  of  6  kW 
provides  a  measure  of  the  average  growth  rate.  Fcr  ex¬ 
ample,  the  peak  growth  rate  is  obtained  at  (o/ckw  ~7  2J 
for  the  fundamental  (upper  intersection),  and  saturation 
occurs  at  kazsal^l84  with  an  efficiency  r;~5.57%.  At 
the  third  harmonic,  peak  growth  is  found  for 
(j/ckw^  27.25,  and  saturation  occurs  at  ku.zsat~172 
with  an  efficiency  77^0.81%.  In  view  of  the  shorter  sat¬ 
uration  length  and  lower  efficiency  at  the  harmonic,  it  is 
clear  that  the  third-harmonic  growth  rate  is  a  substantial 
fraction  of  that  at  the  fundamental.  A  similar  con¬ 
clusion  is  obtained  at  the  seventh  harmonic  which 
displays  a  peak  growth  rate  at  u>/ck u.  ~65.5  and  a  satu¬ 
ration  length  kwzM^  172  with  0.39%. 

The  effect  of  the  thermal  spread  is  shown  in  Fig.  2  in 
which  we  plot  the  normalized  efficiency  77/77 0  (770  is  the 
efficiency  for  Ap,  =  0)  versus  the  axial  energy  spread 
Ay,/y o-  Observe  that  since  the  initial  distribution  is 
monoenergetic 


,  Ap, 
l+2(^-l)  — 
P  0 


-1-2 


(12) 


where  Yo=  1  +Pom  The  principal  conclusion  to  be 
drawn  from  the  figure  is  that  the  harmonic  emission  is 
far  more  sensitive  to  the  effect  of  the  thermal  spread 
than  the  fundamental.  While  the  efficiency  decreases  by 
an  order  of  magnitude  for  Ay,/y0>2%  at  the  funda¬ 
mental,  a  corresponding  decrease  occurs  for 


TE01  Mode  (a  =  10  cm;  b  =  3  cm) 


FIG.  2.  Thermal  effect  on  the  efficiency. 


Ay. /y O^0. 7%  at  the  third  harmonic.  At  higher  har¬ 
monics  the  decrease  is  even  more  rapid.  In  addition,  as 
found  by  Davidson,6  the  growth  rates  are  also  more  sen¬ 
sitive  to  the  thermal  spread  at  the  harmonics  than  at  the 
fundamental.  While  the  saturation  length  increases  from 
ka  zsat  ~  199  to  k„zsal  i209  at  io/ckw  =7.0  (fundamental) 
as  the  thermal  spread  increases  to  Ay,/y0^  1%,  the  de¬ 
crease  is  more  rapid  at  the  third  harmonic.  For  this 
case,  the  saturation  length  increases  from  **•**.,  -  -0-  to 
kuzM  -  241  at  io/ckw  =  27. 25  as  the  axial  energy  spread 
increases  to  Ay.  /y0  =  0.5%.  At  co/ck^  =65.25  'seventh 
harmonic)  the  saturation  length  increases  from 
kwz»,  ^205  to  ~  257  as  the  thermal  spread  in¬ 

creases  to  Ay.  /y0^0. 27%. 

The  question  of  the  effect  of  the  thermal  spread  on 
harmonic  generation  has  been  addressed  in  an  analytic 
model  of  the  low-gain,  single-particle  regime  by  Coisson 
and  de  Martini.7  The  results  of  the  analytic  model  also 
showed  that  the  sensitivity  of  the  interaction  to  the 
beam  thermal  spread  increased  with  the  harmonic  num¬ 
ber.  The  reason  for  this  can  be  understood  by  considera¬ 
tion  of  the  resonance  condition  oj  —  (k  +lkw)u!  -0,  so 
that  the  transition  to  the  thermally  dominated  regime 
occurs  when  ( Rek  +  lkw  )Ad,  =  ( Imk )F,  where  Ac,  is  the 
axial  velocity  spread,  £f,  is  the  bulk  axial  velocity  of  the 
beam,  Rek  is  the  wave  number,  and  Imk  is  the  growth 
rate.  Since  the  growth  rate  typically  decreases  with  in¬ 
creasing  harmonic  number6  while  the  wave  number  in¬ 
creases  (i.e.,  it  is  a  higher-frequency  interaction),  it  is 
clear  that  the  transition  to  the  thermally  domirated  re¬ 
gime  occurs  for  progressively  lower  values  of  Ad,  /u,  as 


s 
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the  harmonic  number  increases. 

The  issue  of  beam  propagation  through  the  wiggler 
can  also  be  addressed  with  this  formulation.  Results  in¬ 
dicate  that  while  the  entire  beam  propagates  through  the 
system  in  the  absence  of  a  thermal  spread,  a  substantial 
fraction  of  the  beam  can  be  lost  to  the  wall  for  even 
moderate  thermal  spreads  (i.e.,  Ay:/y0>l%)  at  this 
beam  energy. 

In  summary,  the  results  indicate  that  substantial 
growth  rates  and  efficiencies  at  the  higher  harmonics  are 
obtainable  at  the  cost  of  more  stringent  requirements  on 
beam  quality.  It  should  also  be  remarked  that  while  the 
difficulties  harmonic  emission  posed  for  the  optical  FEL 


oscillator  at  ACO  were  due  to  spontaneous  (rather  than 
coherent)  emission,  the  higher  beam  currents  obtainable 
o-i  the  forthcoming  super-ACO  storage  ring  could  lead 
to  coherent  harmonic  generation.  Such  problems  are 
due,  in  large  part,  to  the  excellent  beam  quality  (i.e., 
Ayz/y0^0.1%)  obtainable  with  a  storage  ring. 
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ABSTRACT 

The  nonlinear  evolution  of  a  Free-Electron  Laser  (FEL)  amplifier  is  investigated 
for  a  configuration  in  which  an  electron  beam  propagates  through  an  overmoded  rect¬ 
angular  waveguide  in  the  presence  of  a  planar  wiggler  with  parabolically  tapered  pole 
pieces.  The  analysis  is  fully  three-dimensional  and  describes  the  evolution  of  an  arbi¬ 
trary  number  of  resonant  TE  and/or  TM  modes  of  the  rectangular  guide  as  well  as  the 
trajectories  of  an  ensemble  of  electrons.  Numerical  simulations  are  conducted  for  pa¬ 
rameters  consistent  with  the  35-GHz  amplifier  experiment  performed  by  Orzechowski 
and  coworkers,  in  which  the  TE0l,TE7i,  and  TM7l  modes  were  observed.  The  the¬ 
ory  is  found  to  be  in  good  agreement  with  the  experiment.  Surprisingly,  comparison 
with  a  single-mode  analysis  shows  that  the  enhancement  of  the  efficiency  of  the  T E0 , 
mode  obtained  by  means  of  a  tapered  wiggler  is  significantly  greater  (as  well  as  being 
in  substantial  agreement  with  the  experiment)  when  the  TE7l  and  TM71  modes  are 
included  in  the  simulation. 
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I.  INTRODUCTION 


The  Free-Electron  Laser  (FEL)  has  been  shown  to  be  a  high  power  radiation 
source  over  a  broad  spectrum  extending  from  microwave1' 13  through  optical  14-21 
wavelengths.  For  operation  at  relatively  low  beam  energies  (typically  below  about  500 
keV)  and  long  wavelengths,  the  device  is  termed  a  Ubitron  1  and  the  interaction  occurs 
in  the  vicinity  of  the  lowest  order  waveguide  cutoff.  As  a  consequence,  the  system 
can  be  designed  in  such  a  way  that  the  beam  is  resonant  only  with  the  lowest  order 
waveguide  mode,  and  a  single-mode  analysis  is  sufficient  to  describe  many  aspects  of 
the  interaction.  However,  at  higher  energies  and  shorter  wavelengths,  the  interaction  is 
overmoded  in  the  sense  that  the  electron  beam  can  be  resonant  with  several  (perhaps 
many)  modes,  and  the  competition  between  the  modes  has  important  consequences  for 
the  interaction. 

The  motivation  for  the  present  work  is  to  develop  a  multi-mode  nonlinear  theory 
and  simulation  code  for  a  Ubitron/FEL  amplifier.  The  analysis  is  based  on  previ¬ 
ously  described  single-mode  analyses  of  helical  wiggler/axial  guide  field22-25  and  pla¬ 
nar  wiggler26  configurations.  The  development  of  a  multi-mode  analysis  represents  a 
straightforward  generalization  of  the  single-mode  theories,  and  involves  the  calculation 
of  J  •  E  for  each  mode  as  well  as  the  integration  of  electron  trajectories  in  the  aggre¬ 
gate  field  composed  of  the  sum  of  all  the  resonant  modes.  The  particular  configuration 
considered  in  the  present  work  is  that  of  a  planar  wiggler  geometry  in  which  the  elec¬ 
tron  beam  propagates  through  a  rectangular  waveguide.  The  detailed  wiggler  model 
we  employ  includes  the  effect  of  parabolically  shaped  pole  pieces  in  order  to  provide 
for  electron  focussing  in  the  plane  of  the  bulk  wiggler  motion,  1,27  and  we  model  the 
injection  of  the  electron  beam  into  the  wiggler  by  allowing  the  wiggler  amplitude  to 
increase  adiabaticaily  from  zero  to  a  constant  level.  In  addition,  we  consider  the  effect 
of  a  tapered  wiggler  amplitude  on  efficiency  enhancement  in  overmoded  systems.  As  in 
the  case  of  the  single-mode  analysis,26  the  overlap  between  the  electron  beam  and  the 
transverse  mode  structure  of  either  TE  or  TM  modes  is  included  in  a  self-consistent 


2 


way,  and  no  arbitrary  “filling-factor”  is  necessary.  Although  the  problem  of  interest 
is  that  of  an  overmoded  FEL  amplifier  which  requires  a  multi-mode  treatment,  only 
single-frequency  propagation  need  be  considered.  As  a  result,  Maxwell’s  equations  may 
be  averaged  over  a  wave  period  which  results  in  the  elimination  of  the  fast-time-  scale 
phenomena  from  the  formulation. 

The  organization  of  the  paper  is  as  follows.  The  general  formulation  is  described 
in  Sec.  II,  and  allows  for  the  inclusion  of  an  arbitrary  number  of  modes  of  TE  and/or 
TM  polarization,  subject  to  the  restriction  that  all  are  propagating  modes  at  the 
same  frequency.  A  direct  application  of  the  multi-mode  analysis  is  to  the  description 
of  a  recent  experiment  by.  and  coworkers,9'13  in  which  the  TE01  ,TE2l,  and  TM2i 
modes  of  a  rectangular  waveguide  were  observed.  Numerical  examples  appropriate  to 
this  experiment  are  discussed  in  Sec.  Ill,  and  good  agreement  with  the  experiment  is 
found.  A  summary  and  discussion  is  given  in  Sec.  IV. 


n.  GENERAL  FORMULATION 


The  configuration  we  employ  is  that  of  an  electron  beam  propagating  through  an 
overmoded  rectangular  waveguide  in  the  presence  of  a  planar  wiggler  field  generated 
by  a  magnet  array  with  parabolically  tapered  pole  pieces.1,27  As  a  result,  the  wiggler 
field  is  assumed  to  be  of  the  form: 


B„,  (x)  =  Bw  {cos  kw  z(sinh[-^]  sinh}-^]^  cosh[-^]  cosh[-^]e„) 


V2 


V2 


V2 


V2 


k*y 
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where  Bw  denotes  the  wiggler  amplitude,  and  kv  (=  2n/\w)  is  the  wiggler  wavenumber. 
The  injection  of  the  beam  into  the  wiggler  is  modeled  by  an  adiabatic  increase  in  the 
wiggler  amplitude  over  Nw  periods.  In  addition,  since  the  enhancement  of  the  efficiency 
by  means  of  a  tapered  wiggler  is  also  studied,  the  wiggler  amplitude  will  be  tapered 
downward  starting  at  some  point  zQ  downstream  from  the  entry  region  in  a  linear 
fashion.  For  this  purpose  we  choose 


|  sin  kw  ze, } 


(1) 


Bw  (z) 


'  Bw  sin* (kwz/ANw)  \0<z<NwXw 
•  Bw  \NW\W  <  z  <  zQ 

.  Bw  [1  +  €wkw  (z  -  2o)l  \z>z0 
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where 
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describes  the  slope  of  the  taper.  Since  the  fringing  fields  associated  with  the  tapered 
wiggler  amplitude  are  neglected,  this  representation  requires  the  slopes  of  the  taper  to 
be  small  (i.e.,  Nm  must  be  large  and  |e„  |  <  1). 

The  boundary  conditions  at  the  waveguide  wall  may  be  satisfied  by  expanding  the 
vector  potential  in  terms  of  the  orthogonal  basis  functions  of  the  vacuum  waveguide. 
Thus,  we  write  the  vector  potential  of  the  radiation  in  the  form 


oo  t 

6A(x,t)  =  SAen(z)e(e1n>(x,i/)  cos  a  (4) 

i,n  =  0 


for  the  TE  modes,  and 


<5A(x,f)=  Yl  ^n(2)[e(^’(i,y)cos  a  ^ 

t,n  =  1 


kln  .  ,tnX^  .  ,nnY 


)sin( — - — )sinae,]  (5) 


for  the  TM  modes,  where  for  frequency  u  and  wavenumber  k(z) 


a  =  S  dz'k(z')  —  ut  . 

O 


(6) 


In  addition,  Y'  Indicates  that  t  and  n  are  not  both  zero,  and 


a  /  i )  /  \  n7T  ,iirXx  .  ,nirYs* 
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are  the  polarization  vectors.  In  this  representation,  the  waveguide  is  assumed  to  be 
centered  at  the  origin  and  bounded  by  —a/2  <  x  <  a/2  and  —6/2  <  y  <  6/2.  As  a 
consequence,  X  =  x  +  a/2,  Y  =  y  +  6/2,  and 


ktn 


(9) 


denotes  the  cutoff  wavevector.  It  is  implicitly  assumed  that  both  6Atn(z)  and  k(z) 
vary  slowly  over  a  wave  period. 

The  multi-mode  treatment  includes  an  arbitrary  number  of  propagating  modes  of 
TE  and/or  TM  polarization.  The  detailed  equations  which  describe  the  evolution  of 
the  amplitudes  and  wavenumbers  of  these  modes  are  identical  to  those  derived  in  the 
single-mode  analysis,26  and  we  merely  restate  the  results  here.  The  equations  which 
govern  the  evolution  of  the  T Etn  mode  are 


and 


d2  ,w2  2 


*?.)««<.  =8^F,.<^e  ”’-v 


(10) 


2ki''-(kl/26atn)  =  -8  ^Ftn 


rsina  (1) 
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(11) 


where  6atn  =  c6Atn/mc 2  =  4nnbe3 /m  (where  nb  is  the  bulk  density  of  the  beam), 

v  is  the  instantaneous  electron  velocity,  and  Ftn  =1/2  when  either  t  =  0  or  n  =  0  and 
unity  otherwise.  For  the  TMfn  mode  we  obtain  a  similar  result 
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where  we  note  that  there  is  no  nontrivial  TM  mode  solution  when  both  i  —  0  and 
n  --  0. 


Equations  (10)  -  (13)  are  equivalent  to  a  calculation  of  J  •  5Efn  for  each  mode. 
The  averaging  operator  ((...))  is  defined  over  the  initial  conditions  of  the  beam,  and 
includes  the  effect  of  an  initial  momentum  spread  by  means  of  the  distribution  function 


^o(Po)  =  4exp[-(plo -p„)2/Ap*]  6[p20 -p\0 -p]0)H[pt0)  (14) 

where  p0  and  A p,  describe  the  initial  bulk  momentum  and  momentum  spread,  H(x) 
is  the  Heaviside  function;  and  the  normalization  constant  is 


A  =  (?r /<fp,0exp[-(p,0  -  poy  / Ap]})~ 1  .  (15) 
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Observe  that  this  distribution  is  monoenereetic.  but  contains  a  pitch  angle  spread  which 
describes  an  axial  energy  spread  given  approximately  by 


At, 

lo 


S  1  -  [1  +  2(t l  -  1) 


1/2 


(16) 


where  %  =  (1  +  p^/m2c2)l/'3 .  As  a  result,  the  averaging  operator  takes  the  form 


((•••)>  =  /  d<t>°  /  dPzo0to  exp[— (plo  -  Po)2/Ap*] 

Zirao  o  o 

if  «/2  6/2 

x  /  <ty>o<T||(0o)  /  dzD  /  dy0o±  (i0,y0 )(••■)  (17) 

-ir  -a/3  —  b/2 

where  r/>0(=  ~cji0)  is  the  initial  ponderomotive  phase,  </>„  =  tan ~l  [pyo/pzo),  0,o  = 
uIO/c,  and  <7||  rip0}  and  CTj.  ( x0,y0 )  describe  the  initial  beam  distributions  in  phase  and 
cross  section. 

The  phase  variation  of  each  mode  can  be  analyzed  by  the  addition  of  an  equation 
to  integrate  the  relative  phase. 


*{z)  =  ]dZ'(k(z')-k0)  ,  (18) 

0 

where  ka  =  (wJ/ca  -  k3(n)l/2  is  the  wavenumber  of  the  vacuum  guide.  Since  the 
departure  of  lfc(z)  from  the  vacuum  wavenumber  describes  the  effect  of  the  wave-particle 
interaction,  $(z)  represents  a  measure  of  the  dielectric  effect  of  the  FEL  interaction. 
Thus,  we  integrate  the  additional  equation 


8 


(19) 
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4-*  =  k~ 


ka 


for  each  TE  and  TM  mode. 

Each  mode  will  interact  resonantly  with  the  electrons  and  be  coupled  through  the 
electron  motion  in  the  combined  wiggler  and  bulk  radiation  fields.  Thus  in  order  to 
complete  the  formulation,  the  electron  orbit  equations  must  also  be  specified.  Since  we 
describe  an  amplifier  model,  we  choose  to  integrate  in  z,  and  write  the  Lorentz  force 
equations  in  the  form 


v. 


d_ 

dz 


p  =  -e6E  -  -v  •  (B„  +  6B) 


where  B„  is  given  by  Eq.  (1)  and  the  radiation  fields  are  given  by 


(20) 
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and 


<5B  =  Vx  £  6Aen  .  (22) 

alt  mode * 

Finally,  the  electron  coordinates  obey  the  equations 
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(23) 
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which  describes  the  evolution  of  ponderomotive  phase 


in.  NUMERICAL  ANALYSIS 


The  dynamical  equations  for  the  particles  and  fields  described  in  Sec  II  are  now 
solved  for  an  overmoded  amplifier  configuration  in  which  several  modes  may  be  in 
resonance  with  the  beam  at  a  fixed  frequency  u.  The  numerical  problem  involves  the 
solution  of  a  set  of  6 NT  +  4JVW  ordinary  differential  equations  (wherejVr  is  the  number 
of  particles,  and  NM  is  the  number  of  modes)  as  an  initial  value  problem.  Observe  that 
equations  for  the  amplitude,  growth  rate,  wavenumber,  and  phase  are  integrated  for 
each  mode.  The  integration  is  accomplished  by  means  of  a  fourth  order  Runge-Kutta- 
Gill  technique,  and  the  particle  average  described  in  Eq.  (17)  is  performed  by  and  N"* 
order  Gaussian  quadrature  in  each  of  the  initial  variables.  The  initial  conditions  on  the 
fields  are  chosen  to  model  the  injection  of  an  arbitrary  power  level  of  each  mode,  and 
the  initial  wavenumbers  correspond  to  the  vacuum  state  (i.e.,  k{z  =  0)  =  ka).  Further, 
the  initial  value  of  the  relative  phase  of  each  mode  is  zero,  and  both  the  wiggler  field 
and  growth  rate  are  initially  zero.  The  initial  state  of  the  electron  beam  is  chosen  to 
model  the  injection  of  a  continuous,  axisymmetric  electron  beam  with  a  uniform  cross 
section  so  that  o\\  =  1  for  —  n  <  0O  <  n,  and  a±  —  1  for  r„  <  Rb.  A  more  detailed 
description  of  the  procedure  is  to  be  found  in  refs.  22  and  26. 

The  particular  example  we  consider  is  that  of  a  35-GHz  amplifier  employing  an 
electron  beam  with  an  energy  of  3.5  MeV,  a  current  of  850A,  and  an  initial  radius  of 
1.0  cm  propagating  through  a  waveguide  characterized  by  a  =  9.8 cm  and  b  =  2.9 cm. 
Wave- particle  resonance  is  obtained  in  the  vicinity  of  35-GHz  for  a  wiggler  field  of 
3.72  kG  amplitude  and  9.8  cm  period, and  beam  injection  is  accomplished  over  an  en¬ 
try  length  of  five  wiggler  periods.  For  this  choice  of  parameters  three  wave  modes 
are  resonant;  specifically,  the  TE0i ,  TE2i,  and  TM2l  modes.  The  multi-mode  results 
described  herein  will  be  compared  with  a  previous  single-mode  treatment.  In  addi¬ 
tion,  the  parameters  correspond  with  an  experiment  conducted  by  Orzechowski  and 
coworkers,913  and  a  comparison  of  the  simulation  with  the  experiment  will  be  given  in 
Sec.  IV. 
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The  detailed  evolution  of  the  total  wave  power  as  a  function  of  axial  position  is 
shown  in  Fig.  1  for  an  axial  energy  spread  of  A 7*/7o  =  1.5%  on  the  beam  and  the 
injection  of  a  signal  at  w/cfcB  =  11.3(34.6Gffz)  composed  of  the  TE01  mode  at  50 
kW,  the  TE2 1  mode  at  500  W,  and  the  TM21  mode  at  100  W.  As  shown  in  the  figure, 
saturation  of  the  total  signal  occurs  at  kw  z  =  96(1.5m)  at  a  power  level  of  201  MW 
for  a  total  efficiency  of  rj  =£  6.87%.  It  is  also  evident  that  although  the  TE0 1  mode 
was  the  overwhelming  dominant  mode  upon  injection,  it  comprises  only  about  60%  of 
the  signal  at  saturation.  The  remaining  power  is  composed  primarily  of  the  TE21  (at 
37%  of  the  signal)  with  only  a  relatively  small  contribution  of  the  TM0l  mode.  As 
discussed  in  ref.  26,  the  reason  for  this  is  that  at  this  frequency  the  growth  rate  of 
the  TE21  mode  exceeds  that  of  the  TE01  mode,  and  compensates  for  the  lower  initial 
power  level.  Due  to  the  polarization  of  the  TM2l  mode,  the  growth  rate  and  efficiency 
are  smaller  than  for  the  TE  modes,  and  the  TM2[  mode  never  accounts  for  more  than 
about  7  MW.  The  rapid  oscillation  shown  in  the  figure  has  a  period  of  approximately 
Xw  /2  and  occurs  because  the  evolution  of  J  •  E  for  a  planar  wiggler  exhibits  both  a  slow 
variation  corresponding  to  the  ponderomotive  phase  and  a  rapid  oscillation  at  A„,  /'2.26 
Observe  that  the  single-mode  analysis  showed  a  saturated  power  of  162  MW  for  the 
TE0l  mode,  126  MW  for  the  TE21  mode,  and  25  MW  for  the  TM2l  mode.  Thus, 
while  the  total  power  of  the  signal  in  the  multi-mode  analysis  somewhat  exceeds  that 
shown  in  the  single-mode  cases  (for  the  TE  modes),  the  power  levels  of  the  individual 
modes  are  lower. 

The  phase  variation  of  each  of  these  modes  is  shown  in  Fig.  2  as  a  function  of 
axial  position,  where  the  arrow  denotes  the  point  at  which  the  total  power  saturates. 
Of  these  modes,  the  TE0i  behaves  in  qualitatively  the  same  way  as  in  the  single-mode 
case.  Specifically,  the  bulk  phase  at  this  frequency  (apart  from  the  rapid  oscillation  at 
Xw  /2)  increases  monotonically  with  axial  position  through,  and  beyond,  the  saturation 
point.  In  contrast,  the  relative  phases  of  both  the  TE2i  and  TM2\  modes  are  decidedly 
not  monotonic  and  exhibit  a  decrease  with  axial  position  starting  at  a  point  somewhat 
beyond  saturation.  This  is  a  multi-mode  effect  since  the  relative  phases  of  the  TE2l 
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and  TM31  modes  also  exhibit  a  monotonic  increase  with  axial  position  at  this  frequency 
in  the  single-mode  analysis.  Finally,  we  observe  that  the  curves  of  relative  phase  for  the 
TE31  and  TM3l  modes  are  almost  identical.  The  reason  for  this  is  that  the  dispersion 
curves  for  the  TEln  and  TMtn  modes  are  degenerate  in  a  rectangular  waveguide. 

The  effect  of  an  initial  momentum  spread  on  the  saturation  efficiency  of  the  total 
signal  and  the  TE01  and  TE3l  mode  components  is  shown  in  Fig.  3.  The  TM21  mode 
is  excluded  from  the  figure  because  it  composes  such  a  small  fraction  of  the  signal.  As 
shown  in  the  figure,  the  saturation  efficiency  is  relatively  insensitive  to  the  axial  energy 
spread  over  the  range  A7,/70  <2.5%,  and  decreases  from  rj  ~  8.6%  at  A  7,  ~  0  to 
rj  ~  5.9%  at  A 7,/7„  zz  2.5%.  The  reason  for  this  is  that  the  coupling  coefficient  (and, 
hence,  the  growth  rate)  depends  upon  the  product  of  the  wiggler  amplitude  and  period. 
Since  this  product  is  large  for  the  present  choice  of  parameters,  the  growth  rate  is  large 
and  the  interaction  can  accept  a  relatively  large  axial  energy  spread  without  suffering 
a  severe  degradation. 

The  saturation  efficiency  is  known  to  scale  as  the  cube  root  of  the  beam  current  at 
the  frequency  of  peak  growth  from  the  idealized  one-dimensional  theory  of  the  high-g-un 
Compton  (i.e.,  the  strong-pump)  regime,  and  this  type  of  scaling  law  was  also  found 
from  the  three-dimensional  single-mode  simulation  of  this  configuration.  The  scaling 
of  the  total  power  as  a  function  of  beam  current  for  the  multi-mode  analysis  is  shown 
in  Fig.  4  for  A  7,  =  0  and  A7,/70  =  1%,  the  efficiency  is  found  to  scale  approximately 

rl/3 

as  17  ~  /„  . 

Turning  to  the  question  of  the  enhancement  of  the  efficiency  by  means  of  a  tapered 
wiggler,  we  plot  the  evolution  of  the  power  with  axial  position  in  Fig.  5  for  parameters 
corresponding  to  those  shown  in  Fig.  1.  The  optimal  start-taper  point  for  this  case 
is  kmz„  =  86,  and  we  choose  a  slope  of  c„  =  -0.007  which  was  also  studied  for  the 
single-mode  analysis.  The  central  conclusion  to  be  drawn  from  the  figure  is  that  it  is 
possible  to  selectively  enhance  the  T E0l  mode.  The  uniform  wiggler  interaction  for  this 
example  yields  a  total  efficiency  of  6.8%,  of  which  the  TE0 1  mode  comprises  only  60% 
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of  the  signal.  By  contrast,  careful  choice  of  both  the  start-taper  point  and  the  slope  of 
the  taper  show  that  the  efficiency  can  be  enhanced  to  rjmax  ~  41.29%  (if  the  wiggler 
field  is  tapered  to  zero)  with  99%  of  the  power  in  the  TE0i  mode.  Both  the  TE2 1 
and  TM,!  (not  shown  in  the  figure)  ultimately  decay  to  extremely  low  intensities.  One 
surprising  result  of  the  present  multi-mode  analysis  is  that  the  maximum  efficiency  to  be 
obtained  by  tapered  wiggler  fields  is  enhanced  relative  to  the  single-mode  analysis.  By 
comparison,  the  single-mode  analysis  for  these  parameters  yields  a  maximum  efficiency 
of  rjmax  ~  34%  which  is  substantially  lower  than  the  41.29%  found  in  the  multi-mode 
simulation.  The  phase  variation  of  the  TE01  mode  for  this  example  is  shown  in  Fig. 
6,  and  exhibits  the  same  qualitative  variation  as  in  the  single-mode  analysis.  Another 
characteristic  of  the  tapered  wiggler  interaction  observed  in  the  single-mode  treatment 
is  that  the  overall  efficiency  appears  to  be  relatively  insensitive  to  the  axial  energy 
spread.  As  shown  in  Fig.  7,  in  which  we  plot  the  maximum  obtainable  efficiency  versus 
A  7,  / 70 ,  this  is  also  found  to  be  the  case  in  the  multi- mode  analysis.  As  shown  in  the 
figure,  the  maximum  efficiency  decreases  from  43.6%  at  7,  =  0  to  as  much  as  39.4%  at 
A' uho  =  2%.  This  is  a  much  lower  proportional  sensitivity  to  the  axial  energy  spread 
than  is  illustrated  in  Fig.  3  for  the  uniform  wiggler  case. 

Finally,  we  address  the  question  of  the  sensitivity  of  the  tapered  wiggler  interaction 
to  fluctuations  in  the  bulk  energy  of  the  beam.  The  reason  for  concern  with  this 
issue  is  that  the  tapered  wiggler  interaction  is  known  to  be  sensitive  to  the  start-taper 
point.  In  particular,  the  taper  should  begin  at  a  point  shortly  prior  to  saturation 
(for  the  untapered  wiggler)  corresponding  to  the  trapping  of  the  bulk  of  the  beam  in 
the  ponderomotive  potential.  Changes,  or  fluctuations,  in  the  beam  energy  at  fixed 
frequency  are  equivalent  to  variation  of  the  frequency  at  fixed  energy,  and  result  in 
shifts  in  the  growth  rate  and  saturation  point.  For  this  reason,  it  might  be  expected 
that  the  tapered  wiggler  interaction  is  sensitive  to  fluctuations  in  the  bulk  energy  of  the 
beam.  In  order  to  address  this  question,  the  variation  in  the  efficiency  has  been  studied 
as  a  function  of  beam  energy,  and  the  results  are  shown  in  Fig.  8.  For  convenience,  this 
figure  has  been  generated  for  the  limiting  case  of  zero  axial  energy  spread  for  which  the 
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optimal  start-taper  point  is  kw  z0  =  83  at  a  beam  energy  of  3.5  MeV.  Hence,  choosing 
=  -0.007  and  the  aforementioned  start-taper  point,  Fig.  8  describes  the  variation  in 
the  efficiency  with  beam  energy  when  (1)  the  wiggler  is  tapered  to  zero  (A Bw /Ba  =  1), 
and  (2)  when  the  wiggler  is  tapered  to  half  its  ambient  level  (A Bn/B„  =  0.5).  As 
shown  in  the  figure,  there  is  a  sharp  decline  in  the  efficiency  above  approximately  3.55 
MeV.  In  contrast,  there  is  a  more  gradual  decrease  in  the  efficiency  for  energies  down  to 
3.3  MeV,  below  which  the  resonant  interaction  at  ui/ckw  =  11.3  is  lost.  As  a  result,  the 
tapered  wiggler  interaction  will  tolerate  a  bulk  energy  fluctuation  of  the  order  of  8.6  % 
for  these  parameters  without  severe  degradation  in  performance  for  these  parameters. 
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IV.  SUMMARY  AND  DISCUSSION 


In  this  paper,  a  multi-mode  analysis  and  simulation  of  FEL  amplifiers  in  three- 
dimensions  has  been  given  for  a  configuration  in  which  a  relativistic  electron  beam 
propagates  through  an  overmoded  rectangular  waveguide  in  the  presence  of  a  planar 
wiggler  generated  by  means  of  an  array  of  magnets  with  tapered  pole  pieces.  The 
multi-mode  analysis  is  accomplished  by  expansion  of  the  radiation  field  in  terms  of 
the  vacuum  waveguide  modes,  and  an  arbitrary  number  of  propagating  TE  and/or 
TM  modes  is  included  in  the  analysis.  Although  multiple  modes  are  included  in  the 
analysis,  the  problem  of  interest  is  that  of  an  amplifier  and  single  frequency  propagation 
is  considered.  As  a  result,  the  field  equations  are  averaged  over  a  wave  period  in  order 
to  eliminate  the  fast-time  scale  phenomena.  However,  no  average  of  the  orbit  equations 
was  performed,  and  the  electron  dynamics  were  treated  by  means  of  the  fully  three- 
dimensional  Lorentz  force  equations.  As  a  result,  the  effects  of  the  adiabatic  injection 
process,  bulk  wiggler  motion,  Betatron  oscillations,  velocity  shear,  beam  focussing 
due  to  the  wiggler  gradients,  and  phase  trapping  of  the  beam  in  the  ponderomotive 
potential  formed  by  the  beating  of  the  wiggler  and  radiation  fields,  are  all  included  in 
a  self-consistent  manner. 

The  numerical  example  describes  a  35-GHz  amplifier  which  employs  a  3.5 
MeV/850A  electron  beam  with  a  1.0  cm  initial  radius  propagating  through  a  rect¬ 
angular  waveguide  with  dimensions  a  =  9.8cm  and  b  =  2.9cm  in  the  presence  of  a 
wiggler  field  with  a  3.72  kG  amplitude  and  9.8cm  period.  Three  distinct  wave  modes 
are  found  to  be  resonant;  specifically,  the  TE0i,TE2l,  and  TM2i  modes.  The  simula¬ 
tion  is  carried  out  under  the  assumption  that  the  injected  signal  consists  primarily  of 
the  TEoi  mode  at  a  50kW  power  level,  the  TE2l  mode  at  500W,  and  the  TM21  mode 
at  100 W.  Results  indicate  that  although  the  TE2l  mode  was  at  a  relatively  low  initial 
power  level,  it  comprises  upwards  of  3796  of  the  saturated  signal.  The  coupling  between 
the  beam  and  the  TM2l  mode  was  weaker  than  for  the  TE  modes,  and  never  accounted 
for  more  than  a  few  percent  of  the  total  signal.  Comparison  with  a  previous  single-mode 
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analysis38  indicates  that  the  efficiency  of  the  total  signal  is  somewhat  higher  than  that 
found  for  single  modes  in  the  case  of  a  uniform  wiggler.  A  more  dramatic  difference 
between  the  multi-mode  and  single-mode  treatments  is  found  for  a  tapered  wiggler 
interaction.  In  this  case,  it  is  found  that  the  selective  enhancement  of  the  TE0i  mode 
is  possible  and,  indeed,  has  been  experimentally  observed.13  However,  the  power  levels 
to  be  obtained  in  the  TE0i  mode  through  the  multi- mode  tapered  wiggler  interaction 
were  found  to  be  substantially  higher  than  found  in  the  single-mode  simulation.  This 
constitutes  an  important  question  for  future  study. 

The  configuration  and  parameters  described  in  this  paper  nominally  correspond  to 
the  experiment  performed  by  Orzechowski  and  coworkers.9'13  The  principal  differences 
between  the  analytical  configuration  and  the  experiment  are  that  in  the  experiment 
(l)  the  beam  was  injected  into  the  wiggler  through  an  entry  taper  region  one  wiggler 
period  long,  and  (2)  a  quadrupole  field  was  used  to  provide  additional  electron  focussing 
instead  of  parabolically  tapered  pole  pieces.  Since  the  fringing  fields  associated  with 
the  wiggler  field  in  the  entry  taper  region  are  not  included  in  the  analytical  model,  it 
would  bt  invalid  to  apply  the  analysis  for  Nw  =  1.  However,  the  choice  of  Nw  =5 
is  made  as  a  compromise  and  gives  good  agreement  with  experiment  subject  to  the 
additional  assumption  of  an  axial  energy  spread  of  A -7,  /~f0  =  1.5%.  This  is  within 
an  upper  bound  of  2%  on  the  axial  energy  spread  established  by  means  of  an  electron 
spectrometer  measurement.38  The  experimental  measurement  for  a  uniform  wiggler 
interaction  resulted  in  a  saturated  power  level  of  180  MW  over  a  length  of  1.3m.  As 
shown  by  Fig.  1,  the  simulation  gives  a  peak  power  of  204MW  which,  if  we  average  over 
the  fast  Aw/2  oscillation,  is  reduced  to  185  MW.  Given  the  experimental  uncertainties 
in  high  power  measurements,  the  latter  figure  is  more  relevant  for  comparison  and 
is  in  substantial  agreement  with  the  experiment.  The  saturation  length  found  from 
simulation  (that  is,  the  length  of  the  uniform  wiggler  region  plus  one  wiggler  period  to 
account  for  the  entry  taper  region)  is  1.1m,  which  is  also  in  good  agreement  with  the 
experiment.  Note  that  rapid  oscillation  in  the  power  and  relative  phase  at  a  period  of 
A„  / 2  is  likely  to  introduce  a  10-20%  uncertainty  in  the  measurement  of  these  quantities. 
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A  comparison  cam  also  be  made  with  the  tapered  wiggler  experiment13  in  which  the 
wiggler  field  was  decreased  by  55%  (A Bm/Bw  =  0.55)  over  a  length  of  1.1m  (i.e., 
=  —0.0078)  and  the  efficiency  was  found  to  increase  to  34%  for  a  total  power  of 
1  GW.  The  evolution  of  the  total  signal  power,  and  that  of  the  TE  modes,  is  shown 
in  Fig.  9  for  parameters  consistent  with  the  tapered  wiggler  experiment  (the  optimal 
start-taper  point  found  in  simulation  was  kwz0  =  86).  As  shown  in  the  figure,  the 
maximum  efficiency  obtained  by  tapering  the  wiggler  field  to  zero  is  approximately 
40.6%,  of  which  more  than  95%  of  the  power  is  contained  in  the  TE0l  mode.  However, 
over  a  length  of  only  1.1m  beyond  the  start-taper  point  (i.e.,  kvz  —  kwza  =  70.5)  the 
efficiency  is  34%,  of  which  approximately  90%  of  the  power  is  in  the  TE0l  mode.  The 
evolution  of  the  relative  phase  for  this  case  is  shown  in  Fig.  10,  in  which  the  relative 
phase  saturates  at  a  value  in  the  neighborhood  of  120°  downstream  from  the  start- 
taper  point.  This  is  in  good  agreement  with  reported  measurements  of  the  evolution  of 
the  relative  phase  in  the  tapered  wiggler  experiment.29  Thus,  within  the  uncertainties, 
imposed  by  the  choices  of  Nw  and  A 'iz >  the  nonlinear  analysis  is  found  to  be  in  good 
agreement  with  the  experimental  measurements  for  both  uniform  and  tapered  wiggler 
interactions. 
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FIGURE  CAPTIONS 


Fig.  1  The  evolution  of  the  wave  power  (both  total  and  TE0l  mode)  with  axial  position. 

Fig.  2  Plots  of  the  evolution  of  the  relative  phase  versus  axial  position  for  (a)  the  TE0l 
mode,  (b)  the  TE7 1  mode,  and  (c)  the  TM21  mode. 

Fig.  3  Variations  of  the  saturation  efficiencies  of  the  total  signal  and  the  T E  modes  versus 
axia.  energy  spread. 

Fig.  4  Graph  showing  the  scaling  of  the  efficiency  of  the  total  signal  with  beam  current 
for  A'if./x  =  0  and  1%. 

Fig.  5  Plot  showing  the  evolution  of  the  total  signal  and  the  TE  modes  for  a  tapered 
wiggler  interaction  characterized  by  ew  =  -0.007  and  kw  za  =  86. 

Fig.  6  Graph  of  the  evolution  of  the  relative  phase  of  the  TE0l  mode  during  the  course 
of  the  tapered  wiggler  interaction. 

Fig.  7  Illustration  of  the  effect  of  an  axial  energy  spread  on  the  tapered  wiggler  interac¬ 
tion.  Observe  that  each  point  corresponds  to  the  optimal  start-taper  point  of  the 
associated  A-/,. 

Fig.  8  The  variation  in  the  efficiency  of  the  tapered  wiggler  interaction  with  fluctuations 
in  the  bulk  energy  of  the  beam. 

Fig.  9  Plot  of  the  evolution  of  the  total  signal  and  the  TE  modes  for  a  tapered  wiggler 
characterized  by  ew  =  -0.0078  and  kK  z0  =  86. 

Fig.  10  Graph  of  the  evolution  of  the  relative  phase  for  a  tapered  wiggler  interaction 
characterized  by  (w  =  -0.0078  and  kw  zQ  —  86. 
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Efficiency  Enhancement  in  Free-Electron  Lasers 
Driven  by  Electromagnetic-Wave  Wigglers 


HENRY  P. 


Abstract—  A  scheme  for  efficiency  enhancement  in  free-electron  la¬ 
sers  which  are  driven  by  electromagnetic-wave  wigglers  is  described 
which  employs  a  tapered  axial  guide  magnetic  Held.  While  tapered 
magnetostatic  wiggler  designs  have  been  proven  to  be  an  effective  means 
of  enhancing  the  efficiency  of  free-electron  lasers,  practical  difficulties 
occur  in  the  design  of  tapered  electromagnetic-wave  wigglers.  For  this 
reason,  it  is  shown  that  a  tapered  axial  guide  Held  can  be  used  in  con¬ 
cert  with  an  electromagnetic-wave  wiggler  for  the  purpose  of  efficiency 
enhancement. 

AN  extensive  body  of  literature  currently  exists  dem¬ 
onstrating  the  application  of  the  free-electron  laser  as 
a  source  of  radiation  over  an  extensive  range  of  wave¬ 
lengths  ranging  from  the  microwave  to  the  optical  spectra 
[  1  ]-[8] .  The  physical  mechanism  depends  upon  the  prop¬ 
agation  of  an  electron  beam  through  a  periodically  rippled 
magnetic  field  referred  to  as  the  wiggler  or  undulator.  Co¬ 
herent  radiation  results  from  the  undulatory  motion  of  the 
electron  beam  which  permits  a  wave-panicle  coupling  to 
the  output  radiation  by  means  of  the  ponderomotive  po¬ 
tential  formed  by  the  beating  of  the  wiggler  and  radiation 
fields.  The  wiggler  field  itself  may  be  either  magnetostatic 
or  electromagnetic  in  nature.  Although  the  bulk  of  exper¬ 
iments  to  date  have  relied  on  magnetostatic  wigglers,  the 
fundamental  principle  has  also  been  demonstrated  using  a 
large-amplitude  electromagnetic  wave  to  induce  the  req¬ 
uisite  undulatory  motion  in  the  electron  beam  [9J.  The 
basic  difference  between  the  two  classes  of  wiggler  lies  in 
the  frequency  of  the  output  radiation,  which  depends  upon 
both  the  wiggler  period  and  the  electron  beam  energy.  In 
the  case  of  a  magnetostatic  wiggler.  the  wavelength  of  the 
output  radiation  scales  as  X  =  X„/2y;  where  \w  denotes 
the  wiggler  period  and  yh  is  the  relativistic  factor  for  the 
beam.  In  contrast,  the  wavelength  of  the  output  radiation 
for  an  electromagnetic-wave  wiggler  scales  as  X  =  X„/ 
4y;.  As  a  result,  for  fixed  wiggler  periods  and  electron 
beam  energies,  the  electromagnetic-wave  wiggler  will 
produce  shorter  output  wavelengths. 

The  question  addressed  in  this  papef  is  a  method  of  en¬ 
hancing  the  efficiency  of  a  free-electron  laser  driven  by 
an  electromagnetic -wave  wiggler.  The  largest  efficiencies 
observed  with  uniform  magnetostatic  wigglers  are  in  the 
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neighborhood  of  13  percent  [10],  [11).  The  term  “uni¬ 
form”  in  this  sense  refers  to  a  wiggler  with  a  constant 
amplitude  and  period.  Although  such  efficiencies  have  not 
yet  been  observed  with  electromagnetic-wave  wigglers, 
the  number  of  such  experiments  is  small  and  there  is  no 
reason  to  think  that  such  efficiencies  are  not  possible. 
However,  for  many  applications,  higher  efficiencies  are 
desirable.  The  theory  of  efficiency  enhancement  by  means 
of  tapered  magnetostatic  wigglers  has  been  amply  dis¬ 
cussed  in  the  literature  [12]-[17]  and  may  be  affected  by 
either  a  tapered  amplitude  or  period.  Recently,  experi¬ 
mental  confirmation  of  the  theory  has  demonstrated  that 
efficiencies  in  the  neighborhood  of  35  percent  are  possible 
using  tapered  magnetostatic  wigglers  [18].  As  in  the  case 
of  magnetostatic  wigglers,  efficiency  enhancement  can 
also  be  achieved  in  free-electron  lasers  using  tapered  elec¬ 
tromagnetic-wave  wigglers.  However,  while  tapered 
magnetostatic  wigglers  are  relatively  easy  to  design  and 
construct,  tapered  electromagnetic-wave  wigglers  present 
technical  difficulties  in  design  due  to  problems  in  mode 
control.  In  particular,  the  coupling  coefficient  between 
modes  in  a  tapered  waveguide  depends  upon  the  slope  of 
the  taper.  Difficulties  may  ensue,  therefore,  if  the  slope 
of  the  taper  is  comparable  to  the  coupling  coefficient  de¬ 
scribing  the  free-electron  laser  interaction. 

For  this  reason,  an  alternate  efficiency  enhancement 
scheme  is  described  in  this  paper  which  relies  on  a  con¬ 
figuration  in  which  a  tapered  axial  guide  field  is  used  in 
addition  to  an  electromagnetic-wave  wiggler  This  ap¬ 
proach  has  the  advantage  of  extreme  ease  of  construction 
Configurations  in  which  an  axial  guide  field  is  used  have 
been  analyzed  both  for  magnetostatic  and  electromag¬ 
netic-wave  wigglers  [19]— [21].  The  axial  guide  field  is 
often  employed  for  the  purpose  of  confining  the  electron 
beam  against  the  effects  of  self-fields,  and  enhancements 
in  both  the  gain  and  efficiency  are  found  both  in  theory 
and  experiment  [3]  for  a  sufficiently  strong  uniform  axial 
guide  field  that  the  Larmor  period  is  comparable  to  the 
wiggler  period.  The  essential  point  which  should  be  em¬ 
phasized  here  is  that  a  tapered  axial  guide  field  can  be 
used  to  enhance  the  efficiency  of  a  free-electron  laser 
without  the  necessity  of  a  strong  guide  field. 

Efficiency  enhancement  by  means  of  a  tapered  guide 
field  has  been  analyzed  for  free-electron  lasers  with  mag¬ 
netostatic  wigglers  [17],  [22).  In  order  to  formulate  the 
problem  for  an  electromagnetic-wave  wiggler.  we  assume 
that  the  axial  guide  field  is  uniform  for  z  <  cn  and  displays 
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a  linear  taper  thereafter.  Thus. 


B0(z) 


S0; 

s0[  1  +  K0U  -  C0)]; 


where  B0  is  a  constant  amplitude  and  «0(  *  d  In  B0(z)/dz ) 
represents  the  inverse  scale  length  for  variation  of  the  ax¬ 
ial  field.  The  electromagnetic  wiggler  is  represented  by  a 
vector  potential  of  the  form 

S„  , 

A„(c,  t)  =  -  —  [e,  cos  (k„z  +  «„./) 

+  ev  sin  ( k„z  +  w„0]  (2) 

where  Bn  denotes  the  amplitude  of  the  w  iggler  magnetic 
field  and  ( a>„,  k„ )  are  the  frequency  and  wave  vector.  Ob¬ 
serve  that  for  positive  u the  Poynting  flux  is  directed 
antiparallel  to  the  c  axis.  The  radiation  field  is  assumed  to 
be  propagating  parallel  to  the  c  axis,  and  is  represented 
by  a  plane  wave 

6A(z.  r)  =  bA(z)  [e,  cos  (kz  -  wt)  -  i,  sin  (kz  -  w/)] 

(3) 


where  oA(z)  denotes  the  slowly  varying  amplitude  of  the 
vector  potential,  w  is  the  frequency,  and  k  is  the  wave 
vector.  The  electron  beam  equilibrium  is  described  by  the 
steady-state  trajectories  of  electrons  in  the  combined  wig¬ 
gler  and  axial  fluid  fields  [19J.  [20J;  hence,  the  equilib- 
r.'im  velocity  is 

!■„  =  r„  [e,  cos  (k„z  +  u >,j) 

*  e,  sin  ( k„ z  +  t ) ]  +  r  e.  (4) 
where  r  is  the  constant  axial  velocity. 


n„  ( u„  +  A'„  v  ) 

k»  [Ho  -  70(  +  kn  v  )] 

f!,,  „  =  i  eB{)  „  /  me  I .  and 


(5) 


(6) 


We  now  determine  the  response  of  the  electron  beam  to 
the  radiation  field,  and  write  r  =  vQ  +  bv  and  >  =  To  + 
by  to  first  order  in  bA.  The  perturbed  orbit  equations  have 
been  solved  for  the  case  of  a  uniform  axial  guide  field 
|2 1 1  Based  upon  this  analysis,  the  perturbed  orbit  equa¬ 
tions  can  be  shown  to  reduce  to 


t  =  K:  ( sin  —  sin 
dz~ 


(7) 


and 


sin 


*o  To 

4>  —  1 

1  4.  ""  *’■' 

o)„  0 „5<? 

<*> 

1  +  - - 

ckw  c 

+  kn  i’i| 

(*  +  kK)c' 


(9) 


In  (8)  and  (9).  ba  =  ebA/mc 2,  0W  &  t\./v, ,.  -y,  =  (1  — 
i’ij /c~ ) “ 1  /2.  and 


Nl- 


M’i 


u-  +  kuv  II 


i  +  lh 
ckw  c  J 


Ti«Oo 


c~  (w„  +  kwi\)‘ 


o0  -  to(w»  + 


10) 


Equation  (7)  describes  the  trapping  of  electrons  in  the 
ponderomotive  potential  formed  by  the  beating  of  the 
wiggler  and  radiation  fields,  and  is  formally  identical  to 
the  dynamical  equation  derived  for  magnetostatic  wig- 
glers.  In  particular,  it  should  be  noted  that  the  term  in  sin 
M  describes  the  bulk  acceleration  or  deceleration  of 
electrons  due  to  the  tapered  guide  field. 

In  order  to  determine  the  implications  of  this  equation 
for  the  phese  trapping  on  the  efficiency  enhancement  in  a 
tapered  system  we  now  calculate  the  small-signal  gain. 
We  assume  the  system  has  run  to  saturation  at  z  =  Co.  at 
which  point  the  tapered  guide  field  begins,  and  we  define 
the  gain  over  a  length  L  relative  to  this  point  as 

„  _  ba[z  =  Z0  +  L)  -  ba(z  =  Z0) 

°L  =  - _ . - •  <  1  1  ) 


As  a  consequence,  it  may  be  shown  that  [21] 


Gl 


y>b  t’„ 

2 kc:  n5a( Co) 


S:o  *•  L 

:o 


dz  (sin  \p) 


(  12) 


where  <  •  •  ■  >  denotes  a  phase  average  over  an  ensemble 
of  electrons.  Since  the  system  is  assumed  to  have  reached 
saturation  at  c  =  Co.  the  bulk  of  the  electron  beam  has 
been  trapped  in  the  ponderomotive  potential  and  \p  =  M 
Hence, 


G,  » 


yo0~i 

2 c:ki  M(c0) 


r  i  -  ♦  * 

i  +  ^  - 

$ 

ck»  c 

k,  ( t»>,  +  t'n ) 

k(<j)  -  u>„) 


(13) 


where  ^  =  (k  +  k„  )c  -  ( w  —  ui„  )t  is  the  panicle  phase 
with  respect  to  the  ponderomotive  wave. 

K:  -  {k  -  k"]  -  0„ba*.  (8) 

ToT't’n 


The  efficiency  enhancement  Aij(L)  can  be  calculated  by 
taking  the  ratio  of  the  increase  in  the  Poynting  flux  to  the 
power  flux  of  the  electron  beam.  Since  the  Poynting  flux 
increases  by  an  amount  AS  =  (uk /2x)GLbA1(z  =  Co), 
the  efficiency  enhancement  is 
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i{L)  =  /3j(/to£.)  — — — 

t  |  U||  1  ^(oJh,  -t-  kwv n) 

ckw  c  kwv<i(u  -  ww)' 

Observe  that  the  efficiency  enhancement  for  the  case  of 
the  magnetostatic  wiggler  [22]  is  recovered  in  the  limit  as 
ww  -*  0. 

In  order  to  obtain  estimates  of  the  efficiency  enhance¬ 
ments  possible,  we  consider  parameters  consistent  with 
those  discussed  previously  [21],  [22].  Specifically,  we 
choose  yo  =  3.5,  Uw/y0ckw  =  0.05,  ub/y'0 /2ckw  =  0.1, 
and  we  assume  that  X*.  =  1  cm.  We  first  consider  Group 
I  orbits  for  which  fl0  <  -y0(w*  +  kwvtl)  and  we  let 
®o/yockw  —  0.5.  This  corresponds  to  an  axial  magnetic 
field  of  18.7  kG.  In  order  to  obtain  the  frequency  of  the 
electromagnetic  wiggler,  we  must  solve  the  electron  orbit 
equations  ((5)  and  (6)]  together  with  the  dispersion  equa¬ 
tion  for  an  electromagnetic  wave  in  a  uniformly  magne¬ 
tized  plasma  [21].  For  this  choice  of  parameters,  we  ob¬ 
tain  uw/ckw  =  1.007  (i.e.,  30.2  GHz),  and  the  interac¬ 
tion  occurs  at  a  resonant  frequency  oj /ckw  =  44.57  with 
a  wavelength  of  0.22  mm.  The  orbit  parameters  are  (3,,  = 
0.956,  vw/c  —  0.067,  and  $  =  1.018.  As  a  consequence, 
the  efficiency  enhancement  per  unit  length  is  approxi¬ 
mately  A i](L)/L  ~  -0.035  k0.  Thus,  if  we  taper  the 
axial  field  downward,  then  we  obtain  an  efficiency  en¬ 
hancement  of  3.5  percent  if  we  taper  the  field  to  zero  ( i.e. . 
k0L  =  -  1 ).  As  the  resonance  is  approached  more  closely, 
the  efficiency  enhancement  increases.  To  see  this,  wt  ob¬ 
serve  that  if  we  double  the  magnetic  field  (B0  =  37.5  kG 
and  Q0/y0ckw  =  1.0)  while  remaining  on  Group  I  frajec- 
tories,  than  we  find  that  u w/ckw  =  1.010  (30.3  CHi), 
c o/ckw  =  41.66  (  X  =  0.24  mm),  fri  =  0.9529,  u„/c  = 
0. 1019,  and  <l>  =  1.117.  Hence,  At ;(£.)/£,=  -0. 10  *0, 
and  an  efficiency  enhancement  of  10  percent  is  possible  if 
the  axial  field  is  tapered  to  zero.  The  efficiency  enhance¬ 
ment  continues  to  increase  as  the  resonance  at  I)0  ~  yo(u„ 
+  kn  v i )  is  approached  further;  however,  v< i  and  the  res¬ 
onant  interaction  frequency  decrease  rapidly  and  it  be¬ 
comes  increasingly  difficult  to  obtain  short  wavelength 
output  radiation.  Finally,  the  mechanism  is  also  operative 
for  Group  II  trajectories  where  Qn  >.  7o(  t»( ),  but 

at  the  cost  of  prohibitively  high  magnetic  fields. 

While  the  specific  theory  of  efficiency  enhancement  de¬ 
scribed  here  is  directly  applicable  only  to  the  strong-pump 
(or  Compton)  regime  of  operation  in  which  space-charge 
fields  are  negligible,  the  fundamental  mechanism  is  also 
appropriate  to  the  collective  Raman  regime  as  well.  The 
principal  advantage  to  the  tapered-guide  field  approach  to 
efficiency  enhancement  is  the  technical  simplicity  with 
which  a  tapered-guide  field  can  be  controlled.  This  con¬ 
trasts  with  the  difficulties  inherent  in  the  control  of  a  ta¬ 
pered  electromagnetic -wave  wiggler.  It  should  be  re¬ 


marked  that  the  effect  of  the  tapered  guide  field  is 
resonantly  enhanced  when  Q0  ~  +  kwv ,)  which, 

depending  upon  the  specific  experimental  parameters,  may 
require  a  strong  axial  field.  However,  the  efficiency  en¬ 
hancement  mechanism  is  operable  at  all  levels  of  the  axial 
field,  and  does  not  require  strong  guide  fields  near  reso¬ 
nance. 
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